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Preface 


This  thesis  is  concerned  with  investigating  the  feasibility  of 
remote  sensing  of  various  types  of  atoms  and  molecules  by  using  laser 
induced  resonance-fluorescence  techniques  from  an  orbiting  platform. 

A  laser,  if  tuned  to  the  correct  frequency,  will  excite  a  resonance 
line  of  an  atom  or  molecule.  If  the  particle  reradiates,  it  does  so  at 
a  frequency  that  is  characteristic  of  that  particle.  Detection  of  this 
signature  allows  verification  of  the  presence  of  the  particular  species. 

I  hope  that  this  work  will  be  useful  in  evaluating  the  feasibility  of 
laser  remote  sensing  from  space  and  hopefully  my  master's  degree  will 
bring  me  one  step  closer  to  obtaining  my  ultimate  goal  of  flying  in 
space  and  perhaps  performing  this  very  experiment. 

I  would  like  to  thank  my  advisor,  MAJ  Langcfor  sharing  his  extensive 
knowledge  of  this  field  with  me  and  for  being  so  patient  and  pleasant 
during  our  work  together.  I  would  also  like  to  thank  my  wife,  Mary  Ann, 
and  my  children,  Wendy  and  Eric  for  putting  up  with  all  of  the  problems 
that  a  typical  tour  of  AFIT  presents.  Finally,  I  would  like  to  thank 
God  for  the  inner  strength  necessary  to  keep  going  and  never  giving  up 
during  these  past  eighteen  months. 


Glenn  Kweder 


Abstract 


Remote  sensing  from  the  space  shuttle  utilizing  laser-induced 
fluorescence  techniques  was  investigated.  A  feasibility  study  of 
LIDAR  (Light  Detection  and  Ranging)  techniques  to  detect  twelve  gases 
of  scientific  and  military  interest  was  carried  out.  Evaluation  of  a 
European  Space  Agency  LIDAR  hardware  design  was  also  a  task.  The 
main  emphasis  was  on  determination  of  background  noise,  spectral 
transitions,  atmospheric  transmission  and  molecular  absorption 
cross-section  of  the  target  in  order  to  determine  the  minimum  target 
concentration  required  to  be  detectable.  Atomic  sodium,  molecular 
iodine  and  hydrogen  flour ide  were  determined  to  be  promising  candidates 
for  LIDAR  resonance  detection.  Other  molecules  investigated,  such 
as  hydrogen  chloride,  xenon  and  krypton  were  considered  marginally 
detectable  at  best  because  of  their  small  absorption  cross-sections. 

Lack  oi:  sufficient  data  prevented  extended  analysis  of  the  remaining 


molecules . 
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I.  Introduction 


In  recent  years,  LIDAR  (light  detection  and  ranging)  has  become  a 
valuable  tool  in  remote  sensing  of  atoms  and  molecules  in  the  environ¬ 
ment.  With  the  introduction  of  tunable  dye  lasers,  it  has  become  pos¬ 
sible  to  excite  certain  molecular  resonances  of  atoms  or  molecules 
(hereafter  referred  to  as  particles)  and  monitor  their  emission  to  det- 
mine  their  existance,  concentration,  and  spatial  extent  in  the  region 
being  investigated.  In  many  instances,  particular  areas  to  be  probed 
are  virtually  inaccessible  to  researchers  because  of  their  remoteness 
of  location  or  their  location  in  hostile  countries.  In  these  types  of 
scenarios,  LIDAR  monitoring  has  the  potential  to  obtain  the  information. 
Other  techniques  exist  such  as  passive  infrared  spectrometry  that  can 
at  times  collect  the  same  data  using  better  developed  technology.  The 
advantage  that  the  LIDAR  has,  however,  is  that  it  can,  in  certain 
instances*  detect  much  smaller  concentrations  of  a  particular  particle. 
The  laser's  narrow  spectral  linewidth  allows  one  to  excite  and  monitor 
a  very  small  spectral  interval  and  hence  minimize  unwanted  background 
from  reflected  and  atmospheric  emissions.  With  the  advent  of  the  space 
shuttle,  a  new  dimension  in  LIDAR  remote  sensing  is  available.  Virtual¬ 
ly  any  part  of  the  world,  any  country,  friend,  or  foe  can  be  surveyed 
within  minutes  or  a  few  hours  for  certain  types  of  particles  to  obtain 
a  deeper  insight  into  natural  processes  or  a  country's  commercial  or 
military  intentions. 

Several  important  issues  must  be  addressed,  however,  before  this 
space-borne  concept  can  become  reality.  Are  tunable  lasers  available 
that  are  powerful  enough  to  excite  molecules  hundreds  of  kilometers 
away,  and  is  the  fluorescence  return  of  the  molecules  strong  enough  to 
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be  detected  over  the  sometimes  overwhelming  noise  background?  Can  a 
laser  beam  remain  relatively  collimated  at  great  distances  so  that  all 
of  the  laser  photons  are  placed  at  the  point  of  interest?  What  is  the 
atmospheric  transmission  of  laser  and  fluorescence  radiation?  Lastly, 
and  most  importantly,  how  well  do  atoms  and  molecules  of  interest 
absorb  and  emit  radiation? 

The  purpose  of  this  thesis  is  to  perform  a  feasibility  study  on  the 
use  of  fluorescent  LIDAR  techniques  from  the  space  shuttle  to  detect 
several  candidate  atoms  and  molecules  of  scientific  and  military  interest 
(see  Table  I).  A  European  Space  Agency  LIDAR  design  was  also  evaluated. 

Table  I 

Candidate  Molecules  for  Investigation 

Molecule  Reason  for  Interest 

Prove/disprove  theory  of 

!.  Sodium  thunderstorm  updraft 

as  a  source  of  upper 
atmospheric  sodium. 


2. 

HC1 

3. 

HF 

4. 

B2°3 

Possible  by-products 
in  rocket  engine 

5. 

hbo2 

combustion. 

6. 

A12°3 

7. 

»2 

8. 

HT 

•  Possible  by-products 

9. 

HTO 

from  a  nuclear  facility 

10. 

*2 

By-product  isotopes 

II. 

Xe 

from  a  nuclear 

12. 

Kr 

explosion. 

By-product  isotopes 
from  a  nuclear 
explosion. 


II .  Theory 

The  Basic  LIDAR  Equation 

LIDAR  resonance  fluorescence  involves  laser  excitation  of  a  reson¬ 
ance  line  in  an  atom  or  molecule  followed  by  detection  of  photons  created 
by  spontaneous  radiative  decay.  The  wavelength  of  the  fluorescent 
photons  can  indicate  the  exact  species  emitting.  For  a  given  scenario, 
it  is  imperative  that  a  sufficient  number  of  photons  reaching  the  detect¬ 
or  be  distinguished  from  the  background  and  detector  noise.  The  avail¬ 
able  number  of  absorbing  particles  is  given  by 

K  =  aJn'(Z)dZ 

where : 

4 

a  =  area  of  laser  beam  coverage  at  target 

n'(Z)  =  number  density  of  absorbers  at  an  altitude  Z  residing  in 

the  proper  energy  state  for  laser  beam  absorption.  This 
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integral  is  known  as  the  column  density  and  has  units  of  (length  ). 

The  area  (a)  of  laser  beam  coverage  is  governed  by  the  beam  diverg¬ 
ence  (0^)  of  the  laser  (see  equation  35).  In  order  to  minimize  noise, 
the  field  of  view  (FOV)  of  the  LIDAR  receiver  is  matched  to  the  beam 
divergence . 

The  limits  on  the  above  integral  are  governed  by  the  sample  time 
(At)  and  delay  time  (t^),  These  will  b’e  defined  later.  The  number  of 
photons  per  unit  area  incident  at  the  region  of  interest  is  given  by: 
photons  per  unit  area  incident  =  — 

assuming  a  uniform  beam, 
whe  re : 

I  =  total  number  of  photons  incident  at  the  target.  The  total 


number  of  photons  absorbed  is  given  by: 


Total  Number  Absorbed  =  IKo 

where: 

a  =  effective  absorption  cross-section 
This  assumes  an  optically  thick  target,  i.e.  n'aAZ«l  where  AZ  is  the 
path  length  of  absorbers. 

The  number  of  photons  radiated  per  unit  solid  angle  is  equal  to 
the  number  absorbed  divided  by  4ir  steradians.  The  number  of  photons 
actually  reaching  the  receiver  in  the  absence  of  atmospheric  transmission 
losses  is  equal  to  the  number  emitted  into  the  solid  angle  subtended  by 
the  receiver  aperture: 

Photons  reaching  receiver  =  - j 

4irZ 

where 

Z  =  range  between  cloud  and  receiver 

A  =  area  of  receiver  aperture 

The  actual  number  of  photons  detected  is  determined  by  three  other  fac¬ 
tors:  detector  quantum  efficiency  (rj)>  fluorescence  efficienty  (Op)  and 
receiver  throughput  efficiency  (e).  Atmospheric  transmission  losses 
result  from  fluorescence  emission  being  absorbed  and  scattered  by  other 
particles  in  the  atmosphere.  If  the  laser  is  at  a  different  wavelength 
than  the  fluorescence,  then  the  loss  factor  is  given  as  the  product  of 
the  transmitted  laser  beam  losses  ( T )  and  the  received  signal  losses 
(t).  Detector  quantum  efficiency  is  the  number  of  electrons  produced 
per  incoming  photon. 

The  value  for  n  is  usually  10%  -  20%  for  detectors  in  visible  or 
infared.  Fluorescence  efficiency,  or  branching  ratio,  is  the  ratio 
of  the  number  of  transitions  that  occur  from  upper  level  (a)  to  a  lower 
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level  (b)  to  the  number  that  decay  from  level  (a)  to  all  possible  lower 
levels.  Receiver  throughput  efficiency  accounts  for  losses  in  reflec¬ 
tion  off  mirrors  and  transmission  through  filters  and  windows.  With  the 
above  factors  included,  one  can  now  write  the  general  I.IDAR  equation 
relating  the  number  of  signal  photons  detected  per  sample  time  to  the 
various  parameters  of  range,  laser  size,  etc: 

N  =  IKaAenOFTRTt  (1) 

4 V 

This  equation  assumes  that  the  emitted  laser  pulse  and  the  fluorescence 
lifetime  (x)  of  the  target  are  short  compared  to  the  integration  time 
(At)  of  the  detector.  Of  the  molecules  considered  in  this  thesis,  this 
was  not  always  the  case,  thus  a  modified  form  of  the  LIDAR  equation  was 
developed  by  the  author.  The  expression  for  the  modified  LIDAR  equation 
will  be  derived  later. 

It  should  be  noted  that  in  most  cases  the  column  density  (K)  is 
expressed  as  the  product  of  an  average  target  number  density  (nr)  and 
column  length  being  sampled  (AZ) .  The  column  length  being  sampled  can 
be  related  to  the  sample  time  by  the  relationship: 

AZ  =  ^cAt  (2) 

where : 

c  =  speed  of  light 

The  longer  the  detector  samples  the  fluorescence  return,  the  larger  is 
the  column  being  sampled. 

The  particular  region  in  the  atmosphere  where  sampling  occurs  is 
governed  by  the  delay  time  (tj)  between  laser  pulse  and  initilization 
of  detector  sampling  (see  fig  I),  l'or  the  case  where  the  shuttle  LIDAR 
is  viewing  straight  down  at  the  earth,  if  t^  is  increased,  then  the 


region  being  sampled  is  farther  away  from  the  shuttle  and  thus  closer 


to  the  earth.  As  t^  is  shortened,  then  sampling  occurs  closer  and 
closer  to  the  shuttle.  In  short,  the  sample  time  (At)  specifies  the 
column  length  being  sampled  while  the  delay  time  (t^)  fixes  the  distance 
between  receiver  and  sampling  region.  Extending  At,  to  increase  the 
number  of  photons  returned  can  be  done  only  to  a  certain  point.  If  the 
cloud  of  interest  is  highly  localized  then  there  is  a  fixed  AZ  where 
molecules  of  interest  reside  and  one  must  adjust  At  by  equation  (2)  to 
obtain  the  optimum  return.  This  statement  is  always  true  if  the  radia¬ 
tive  decay  time  (t)  is  much  shorter  than  At.  If  it  is  not,  then  it 
may  be  advisable  to  extend  At  farther  to  allow  detection  of  photons  from 
late-decaying  particles  but  short  enough  to  avoid  additional  noise  from 
reflected  laser  light  off  of  the  earth. 

The  total  time  from  the  first  laser  pulse  to  the  last  fluorescence 
sampling  is  known  as  (see  fig  1).  If  a  LIDAR  "sounding"  period  is 
very  long  (t^  ~  10  seconds)  then  look  angle  and  range  can  no  longer 
be  treated  as  constants  and  equation  (1)  must  be  modified  (see  fig  2). 

This  scenario  can  get  very  complicated  thus  this  thesis  will  be  restrict¬ 
ed  to  short 

Atmospheric  Transmittance 

The  coefficient  of  transmission  through  3  slab  (AL)  of  atmosphere 
is  defined  as  (ref  1:10): 

tslm  -  <3) 

where : 

Y  =  S  +  S  +  K  +  K  (4) 

mama 

and 
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Sm  =  Molecular  scattering  (Rayleigh) 

S  =  Aerosol  scattering  (Mio) 

3 

K  =  Molecular  absorption  (resonance) 
m 

Ka  =  Aerosol  absorption 

Rayleigh  scattering  takes  place  between  a  photon  and  an  atom  or 
molecule  whose  size  is  smaller  than  the  wavelength  (X)  of  the  photon. 

It's  cross-section  has  a  X  ^  dependence.  Mie  scattering  is  much  more 
complex  and  involves  scattering  of  a  photon  off  of  a  particle  compar¬ 
able  in  size  with  the  wavelengtn  of  the  light.  The  Y  coefficient  is  a 
function  of  altitude  and  wavelength.  Using  the  tables  for  laser  line 
transmittance  (ref  1:17,18),  one  can  determine  the  total  transmittance 
through  many  slabs  of  atmosphere  by  the  product  of  equation  (3)  for  each 
slab  (ref  1 :  12) ; 

TTOT  *  T  =  exp  (~?  W  (5) 

j  J  J 

These  tables  did  not  list  the  desired  laser  wavelengths  used  in  this 
study,  however  it  is  appropriate  to  interpolate  between  two  wavelength 
values  for  Ka>  Sm  and  Sa  because  they  are  slowly  varying  functions. 

This  is  not  the  case  at  all  for  which  is  a  highly  oscillatory  function 
of  wavelength.  Fortunately  in  this  thesis  the  number  density  of  abosrbers 
at  the  wavelengths  of  interest  are  low  enough  to  allow  to  be  dis¬ 
carded. 

Modified  Forms  of  the  Basic  I.IDAR  F.quation 

For  other  than  vertical  sounding,  the  look  angle  (0)  and  range  (Z) 
will  be  larger  than  zero  degrees  and  Zy  respectively  (see  fig  I).  This 
may  happen  frequently  when  the  shuttle’s  orbital  parth  does  not  take 
the  LIDAR  directly  over  the  event  of  interest.  The  transmission  co- 


7 


efficients  (Tt)  and  (TR)  will  take  the  form: 

T(0,A)  =  oxpeXy.AL.)  SGC  0  (6) 

j  J  J 

or: 

T(0,A)  =  TSeC  0  (70 

1 

Where  T  is  given  by  equation  (5).  The  range  at  angle  0  is  given  by  the 
relationship 

Z  =  Zvsec  0  (8) 

where : 

Zv  =  nadir  range  (see  fig  2) 

For  a  sounding  about  a  region  other  than  nadir,  equation  (1)  becomes: 

N  =  lKaA£n0F(TRTt)  SGC  0  (9) 

-  - 

4irZysec  0 

To  repeat,  it  is  assumed  that  xd  is  short  so  that  0  and  Z  can  be  treated 
as  constants.  Also  0  is  assumed  small  in  order  to  use  a  flat  earth 
atmosphere  geometry. 

Equation  (1)  assumed  that  the  lifetime  of  the  state  (t)  was  very 
short  compared  to  the  sample  time  (At)  of  the  fluorescence  received  by 
the  detector.  This  meant  that  it  was  safe  to  assume  that  all  of  the 
fluorescence  photons  from  the  target  emitted  in  the  narrow  solid  angle 
of  the  detector  are  collected  assuming  that 

td  =  2Z/c  (10) 

If  the  lifetime  (x)  is  longer  or  comparable  to  the  sample  time(At)  then 
the  detector  will  not  collect  all  of  the  emitted  photons  but,  in  many 

*  . 

The  proper  delay  time  assures  that  no  photons  emitted  into  the  detector 


field  of  view  are  missed  and  at  the  same  time  minimizes  background  noise 
by  keeping  At  as  small  as  possible. 


cases  only  a  small  fraction  of  them  in  a  sample  time.  If  the  lifetime 
is  long  enough,  it  is  possible,  in  a  later  sample  time,  to  continue 
receiving  photons  from  particles  excited  many  pulses  ago  and  just 
beginning  to  decay.  The  total  number  of  photons  absorbed  per  pulse, 
ignoring  transmission  losses  for  the  present  time,  is  given  by: 

NABS  =  IKa 

where : 

I  =  number  of  photons  per  laser  pulse 

The  total  number  of  photons  decaying  after  the  first  laser  pulse  in  a 
time  interval  (At)  is  given  as: 

N  =  N  -  N  _A* At 

decay  1  ABS  ABS 

where  the  factor  represents  the  number  of  particles  still 

excited  after  a  time  interval  (At).  (A*)  is  the  sum  of  the  Einstein 

coefficient1  for  spontaneous  radiative  transitions  between  the  excited 

and  ground  states  and  a  deexcitation  collision  frequency.  During  the 

second  pulse  period,  the  number  of  photons  absorbed  will  be  the  same  as 

in  the  first  pulse  period  (assuming  that  the  laser  output  is  the  same 

for  each  pulse).  The  number  of  photons  decaying  during  the  second  sample 

time  (At)  after  the  second  laser  pulse  will  be  equal  to: 

N,  o  ?  (N.no  -  NADOe_A*At)  +  N, 
decay  2  ABS  ABS  decay 

2nd  2nd  Lefto/er 

pulse  pulse  1st  pulse 

The  number  of  molecules  decaying  during  the  second  sample  time  will  be 
from  some  molecules  excited  during  the  second  pulse  plus  some  left  over 
from  the  first  pulse  that  are  just  beginning  to  deexcitc.  This  second 
terms  is  equal  to: 

N '  ( 1  -  e~A*At) 

Where  N1  is  the  number  of  excited  molecules  left  over  from  the  first 


pulse  that  are  still  in  the  excited  state  when  the  second  pulse  is 
occuring.  If  a  laser  has  a  pulse  repetition  frequency  of  f  hertz,  Chen 
the  time  between  pulses  (t1)  is  equal  to  1/f^,  This  is  the  amount  of 
time  elapsed  between  laser  pulses  and  represents  the  amount  of  time 
that  the  first  group  of  excited  particles  has  had  to  decay.  Therefore, 
the  number  of  excited  molecules  left  over  the  from  the  first  pulse  at 


the  beginning  of  the  second  sampling  time  is: 


N'  =  N  =  N  e"*t:  + 

ABS  ABS 


1st 

pulse 


1st 

pulse 


Where  represent  respectively  the  radiative  lifetime  of  the  mole¬ 

cule  and  the  time  between  collisions.  If  all  laser  pulses  are  equal 
then  tor-  all  pulses  are  the  same.  The  total  number  of  decaying 

events  during  the  second  sample  time  <^t)  is  then: 

Ndecay  2  *  NABS  ( 1  5  NABSe  ( 1  J 


„  , .  -A*At.  L  -A*t  * . 

N,  q  =  N.__  (1-e  )  (1+e  ) 

decay  2  ABS 

It  is  easy  to  see  that  during  the  third  sample  time,  the  decaying 
particles  will  come  from  molecules  excited  from  the  third  laser  pulse 

in  addition  to  some  left  over  from  the  second  and  first  pulses: 

x,  v,  /,  -A*Atx  t  XT  -A*t 1  -A*Atx 

N,  -x  =■  N.do  (1-e  )  +  N._0e  (1-e  ) 

decay  ■>  ABS  ABS 


-2A*t’  , .  -A*Atx 
*NABSe  °-e  > 


Simplifying,  we  have: 


„  „  ,  ,  -A*At«  -A*t 1  -2A*t\ 

Ndecay  3  =  NABS  °  *  e  >  0  +  6  +  e  > 


The  number  decaying  in  the  Nth  pulse  is  given  by: 
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M  mci  ~A*AtN  ..  A  -A*t’  -2A*t' 

Ndecay  N  =  NABS  °  “  e  >  (l  +  e 


.  +  e 


-(N- 1 )A*t '  ) 


or: 


M  M  c,  -A* AC,  N_!  -  JA*  t ' 

Ndecay  N  =  NABS  (1"e  >  *  e 

J —0 


The  total  number  of  deexcited  molecules  after  N  sample  times  is  given  by 


N ,  =  N  ,  .  +  N ,  N  ,  , 

decay  decay  i  decay  2  decay  N 


Total 

A  little  thought  will  reveal  that  this  expression  can  be  written  as: 

N 


N- 1  A .  ,  _  i 

l  N  (l-e"A5At)  z  -JA*t ' 
Total  i=0  ABS  J=0  6 


decay 


This  can  be  written  more  explicitly  in  the  following  way: 

N- 1  i 

—  A7C  A  r  . 

N 


■decay*  IK°  £ 


E  -JA*/f 
i=0  J=0e  P 


(11) 


Total 

If  the  product  A*t(  is  much  less  than  1,  we  may  expand  the  exponential 
in  a  Taylor  series  keeping  only  the  first  two  terms: 


-JA*  t ' 


*«• » 


e  =  (e_A'C  )  J  =  <1-A*t')  J  (12) 

and  if  A*At  in  equation  (11)  is  small  compared  to  1,  which  is  usually 
the  case,  we  may  say: 

_ A  •S’  A 

(13) 


l_e~A*At  =  A*At 


With  the  above  approximations,  we  can  simplify  equation  (11)  to  read: 

N-l  i 


N,  =  IKdA*At  _ 
decay  E 


z  ('-A*/£P >' 


Total  1=0  J=0 

The  total  number  of  photons  actually  detected  by  the  LIDAR  receiver 


over  N  sample  times  is  given  as: 


N, 


=  _R  Ikcr  (1-e  )  Tr’f ^AriGOp,  £  £  _JA*/f 

Detect  nT  i=0  J-0®  P 


Tot 


4  TIE2 


The  factor  (n^/n^,)  will  be  discussed  shortly. 
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Finally,  the  most  general  ca3e  considers  look  angle  dependence.  Recall¬ 


ing  equations  (7),  (8)  and  (9)  we  have: 


Total  =  Hr  IKoU-e  MneGL  (T  T  )  sec  Q  N  -JA*/f  (14) 

Detect  7T  FRt  Die  p 

DeCect  lxr  i=0  J=0 


i-0  J=0 


47i22sec20 

v 


If  the  conditions  outlined  by  equations  (12)  and  (13)  hold  true,  we  may 

rewrite  equation  (14)  as: 

N  A  N-I  i 

NTotal  =  UR  IKoA*AtAne0..  (T_T  )S£C  U  £  £  ( 1  -A* / f  )J  (15) 

_  .  .  -  J?  K.  C  p 

Detect  *T  i=0  J=0 

2  2 
4liZ  sec  Q 
v 

Quenching 

The  factor  (rtR/n,j,)  occurs  because  only  a  fraction  of  decaying  mole¬ 
cules  decay  radiatively.  Others  decay  non-radiative  by  through  eollision- 
al  deexcitation  with  other  molecules.  This  process  is  known  as  quench¬ 
ing.  Quenching,  or  non-radiative  deexcitation,  as  it  is  sometimes  called 
can  result  in  a  reduction  in  photon  return  if  the  time  between  collisions 
of  molecules  is  comparable  to,  or  smaller  than  the  radiative  lifetime  of 
the  state.  Using  the  relationship: 

jM<u>2=|kT  (16) 

for  a  gas  in  equilibrium,  one  can  solve  for  the  velocity  of  the  mole¬ 
cules.  The  time  between  collisions  is  equal  to  the  mean  free  path 
between  collisions  divided  by  the  velocity  of  the  molecules: 


Tc  =  7^ 


n  o  <u> 
c  c 


where: 


nc  =*  Total  number  density  of  all  molecules  present. 

—  15  2. 

a  =  10  cm  (collisiotial  cross-section  for  quenching),  (ref  22:4217) 
c 
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The  question  now  is:  in  a  sample  time  (At),  how  many  excited  molecules 
will  decay  through  radiative  processes  relative  to  the  total  number  that 
decay  through  both  radiative  and  collisional  processes?  For  a  total 
number  (n)  of  excited  molecules,  the  decrease  of  n  in  a  time  (dt)  is 
given  by: 


where : 


-dn  =  dn_  +  dn„ 
K  L 


(18) 


dn  =  ndt/t 
c  c 


and: 

dnR  =  ndt/tR  C 19) 

By  integrating  over  the  sample  time,  one  can  obtain  the  total  number  of 
molecules  remaining  in  the  excited  state  after  a  time  t: 

dn  =  n  (— —  +  •— )  dt  =*  andt 

TR  tc 


which  yields:  n  =  nQe  at  (20) 

The  total  number  that  have  undergone  a ’transition  (nT>  in  time  t  is  just 
equal  to  (nQ-n) : 


nT  =  nQ  ( 1-e  a^c) 

The  number  of  molecules  which  undergo  a  radiative  transition  in  a  time 
(dt)  are  given  by  equation  (19).  Using  equation  (20),  equation  (19) 


can  be  written  as: 


(21) 


.  n  e_at  dt 
dnR  =  — - 


By  integrating  this  expression  over  time,  the  total  number  of  radiative 
transitions  in  At  can  be  obtained:  Integrating,  we  have: 

n„  t 


/  -.4  / 


e~at  dt 


which  yields: 


n 

o 

*»  - 

R  aT„ 


1-e 


-aAt 


(22) 


In  a  time  (At),  the  fraction  that  have  undergone  a  radiative  decay  is 
given  by: 

,-~aAt 


n  n  1-e 
R  o 


=  1 


°T  aVb  aTR 


This  can  be  written  more  explicitly  as; 


n 


T 


(23) 


From  this  expression,  it  is  easy  to  see  that  as  collisions  become  domin¬ 
ant  (t  decreasing)  then  n _/n,r  decreases  and  visa  versa.  If  there  are 

C  K  1 

no  appricable  collisional  effects,  than  0  and  n  /n_  1  as  ex- 

tv  u  K  1 

pected.  This  just  says  that  all  transitions  are  radiative  ones. 


Effective  Cross-Section 


The  cross-section  (a)  used  in  this  thesis  is  an  effective  cross- 


section  for  the  absorptance  of  a  doppler-broadened  particle  given  a 
doppler  broadened  laser  line.  The  absorption  line  profile  of  the  part- 
* .  icle  is  given  as : 

4  » 
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o(v) 


2  a  (ln2) 
o 


nh  Av' 


exp 


-4(1h  2)  (V~V0)' 

% 


(24) 


where  (ref  2:65): 


e  f  (For  atomic  particles) 
4t0MeC 


AVy  =  linewidth  at  full  width  half  maximum 
e  =  electron  charge 

f  =  oscillator  strength 


£q  =  permittivity  of  free  space 

M  =  electron  mass 
e 

For  molecules  another  form  of  <jQ  is  used  (ref  3): 

2 

g  A'X  .  (For  molecules) 
cr  _  a  ab 

°  ‘  8a8Tr 


(26) 


where: 

g  =  degeneracy  of  level  (a) 

=  degeneracy  of  level  (b) 

A'  =>  Einstein  coefficient  for  spontaneous  transition  between  levels 
(a)  and  (b) . 

Xab  =  wavelength  of  the  transition 
and  finally: 

Avd  =  2V0  (ln2)  *  (27) 

MC4 

where : 

K  >=  Boltzman’s  constant 

T  =  temperature 

M  »  mass  of  particle 

viq  =  line  center  frequency 

The  laser  lineshape  profile  is  similar  to  equation  (24): 
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L(v>  =  2(ln2)^  exp  ~4(In  2) 


n-iAv. 


Av, 


L  “VL 

Note  that  L(v) ,  unlike  cr(v),  is  a  distribution  function.  This  is  be¬ 
cause  the  magnitude  of  the  laser  energy  (1)  is  explicitly  expressed  in 
equation  (I).  The  laser  linewidth  (AvL)  takes  the  same  form  as 
equation  (27)  but  with  different  (T)  and  (M)  values.  In  this  thesis, 
a  typical  laser  linewidth  is  assumed  and  not  calculated.  The  effective 
cross-section  is  a  convolution  of  the  absorption  and  laser  line  pro¬ 
files  over  all  v: 


7 

°eff  =0=/l 


L( v)  o(\>)  dv 


a  = 


o  41n2 
o _ 

7rAvLAvD 


exp  -41n2  (v-VQ)2  (  -~y  +  -~-)dv 


Av; 


Av„ 


(28) 


The  integral  has  the  form: 


Letting: 


OO 

/ 


2  2 

e‘ay  dy 


jn 

a 


a  =  41n2(- 


Avt 


Avr 


equation  (26)  can  be  written  as: 


0  = 


a04(ln2)*a 


ttAv.  Av  a 
L  D 


This  reduces  to: 


In  many  cases: 


a  2qo( ln2)^ 


tt^AVt  Avn( 


1 


I 


LUVD  Av2  Av2 

L  D 


(29) 


Avl  >  >Avd 
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there  fore : 


Avl  Avq  Avj) 


i 

Plugging  this  into  equation  (29)  we  obtain: 

2a  i  In.  2)  ** 

°  =  -  (30) 

TI  A\>l 

Noise 


Noise  can  enter  the  system  in  a  variety  of  ways.  Earth-reflected 
sun,  moon,  and  laser  light  all  result  in  unwanted  photons  reaching  the 
detector.  These  photons  along  with  detector  and  signal  noise  place  a 
lower  limit  on  the  number  of  signal  photons  (N)  that  must  be  counted 
or  detected  in  a  given  sample  time.  This  then  places  a  lower  limit  on 
the  column  density  (K)  oF  absorbers  required  through  equation  (1).  The 
photons  arriving  at  the  detector  are  generated  from  random  processes, 
therefore  they  obey  poisson  statistics.  The  probability  that  m  photons 
will  arrive  in  time  t  is  given  by: 

a  t)n 

P(m,t)  =  — —  exp  (-It)  (31) 

ml  p 


where : 


Ip  =  average  rate  of  photon  arrival 


From  the  poisson  distribution,  the  average  number  of  photons  Cm) 


arriving  in  a  time  t  is  given  by  Ipt  and  the  variance  in  m  is  given  by 
Ipt.  The  signal  to  noise  ratio  is  given  by: 


S 

N 


~  =  —2—,  =(T  t ) 1  =  m 

°p  (I  &  p 

P 


(32) 


where : 


0p  =  square  root  of  the  variance 
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The  above  expression  considers  noise  only  as  a  result  of  the  random 

arrival  of  signal  photons.  The  value  however  must  include  further 

noise  contributions  from  the  sources  mentioned  above.  When  these  are 

included  we  have:  ^ 

a  =  (  Nf  [  ,+N  r  , f  +N  •[  f  +N  t  ,  f  +N  ,t  ,  f  )  (33) 

up  p  d  Eld  p  Mld  p  A 1  d  p  dld  p' 

where: 


f  =  pulse  repetition  frequency 
=  total  sounding  time 

N  «  average  number  of  signal  photons  counted  per  sample  time  (equation  l) 
N  =  average  number  of  photons  counted  from  earth-reflected  sunlight  per 
sample  time. 

=  average  number  of  photons  counted  from  earth  reflected  moonlight 
per  sample  time. 

Na  =  average  number  of  photons  per  sample  time  counted  from  airglow  in 

the  atmosphere  (sodium  case  only) 

=  average  number  of  thermal  electrons  created  per  sample  time. 

It  is  clear  that  the  product  (f^^)  yields  the  total  number  of 

pulses  or  sample  periods  (At)  taken  in  the  total  sounding  period  (t^) • 

Therefore,  the  product  of  N,Ng,N^,NA  and  with  yield  a  total 

number  of  photons  collected  in  the  sounding  period  (x^).  Dark  current 

results  from  thermally  emitted  electrons  from  the  photocathode  of  a 

photomultiplier  tube.  The  rate  of  electron  emission  per  unit  area  is 

given  by  the  Richardson  equation  (ref  9): 

i  =  120T2  exp(^ (34) 

cm  ' 

where : 

T  =  temperature  (°K) 
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4»  =  work  function  of  the  photocathode  material  (cv).  The  amount  of 
current  is  proportional  to  the  area  of  the  photocathode.  Finally,  we 
have : 


where: 

5  »  gain  of  individual  dynode  in  photomultiplier. 

Of  course,  if  some  detector  other  than  a  photomultiplier  is  used,  the 

factor  (&~— ) 2  is  not  used  in  equation  (36).  Also,  if  the  detector  temp- 

6 

erature  is  lowered  enough,  can  be  neglected.  Lastly,  will  only 
occur  when  sodium  LIDAR  soundings  are  considered. 

Values  for  and  are  derived  from  a  radiance  vs.  wavelength 
chart  (ref  4:42).  The  radiance  (L)  lor  a  body  is  the  number  of  photons 


per  second  per  unit  area  per  steradian  emiLtcd  in  a  wavelength  interval 


Equation  (36)  can  bo  solved  for  the  minimum  number  of  signal  photons 
required  per  sample  time  to  achieve  a  signal  to  noise  ratio  (R) . 

Squaring  both  sides  of  equation  (36)  we  obtain: 

N^f  t  ( 

2  _  pTd  (6-1) 

"  (N+N+N  +N,+N  ,)6 
L  rl  A  a 

Rea  .'ranging  to  a  more  familiar  form  we  obtain: 

2  2 

fPTd  N  +(6=I)  n+  (PiT  (WW  =  0 

This  is  a  quadratic  equation  in  N.  Using  the  quadratic  formula  we 

obtain:  r  .  *|. 

„2.  2  4f  t  ,R  6(Nr+NM+NA+N.) 2 

N.  +Ri__+_J Li - L-M  ■■  A-.-d.  .  (41) 

mltl  (6-1)  1(6-1) _ (5-1) _ „  __ 

2f  Tj 

pM 

The  result  of  equation  (41)  can  be  substituted  in  for  N  in  equation  (1) 
and  a  minimum  column  density,  (K)  can  be  obtained  for  a  given  R. 


Laser  Beam  Divergence 

All  of  the  LIDAR  equations  discussed  so  far  have  assumed  that  the 
area  of  the  laser  beam  at  the  target  is  smaller  than,  or  equal  to  the 
area  of  the  cloud  being  investigated.  If  the  laser  beam  were  bigger, 
then  many  laser  photons  would  never  find  a  target  particle  and  the 
return  signal  would  be  smaller.  This  could  be  disasterous  in  many 
sitatuions  where  LIDAR  detection  of  various  molecules  is  marginal  at 
best.  It  is  imperative,  therefore  to  assure  that  the  beam  is  always 
smaller  than  the  target.  This  can  be  accomplished  ironically,  by 
initially  expanding  the  beam  at  the  laser  transmitter.  The  full  angle 
beam  divergence  is  given  from  Fraunhofer  diffraction' theory  to  be: 
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(42) 


«  beam  divergence  in  radians 
D  =  diameter  of  laser  beam  exit  aperture 

For  small  beam  divergences  we  may  approximate  the  diameter  (d)  of  Che 
laser  beam  spot  at  a  distance  (Z)  from  the  shuttle  LI DAK  platform  by 
the  relationship: 

d  =  0DZ  (43) 

Using  equation  (42)  we  may  rewrite  equation  (43)  to  directly  relate 
aperture  diameter,  range  and  beam  spot  size: 


D 


1.22XZ 

d 


(44) 


Range,  laser  wavelength  and  target  extent  will  have  to  be  evaluated 
for  each  scenario  considered  to  determine  if  the  laser  on  board  the 


shuttle  needs  an  expanding  telescope  or  not 


III.  LIDAR  Conceptual  Design 


The  LIDAR  hardware  consists  of  two  main  parts.  The  laser  trans¬ 
mitter,  and  receiver  (fig  3).  The  laser  transmitter  in  most  cases  is  a 
tunable  dye  laser  with  expanding  optics  if  needed.  The  receiver  could 
be  a  Cassegrain  telescope  with  folding  mirrors  to  minimize  the  extent 
of  the  optical  train  (ref  5:191).  In  order  to  obtain  a  very  narrow 
spectral  window,  (on  the  order  of  Angstroms  or  fractions  of  an  angstrom), 
a  Fabry-Perot  or  several  Fabry-Perot  etalons  in  series  are  used.  In 
order  to  assure  proper  functioning  of  the  etalons,  a  collimating  lens 
is  used  ahead  of  the  etalons.  The  transmitted  signal  through  the  etalon 
is  then  sent  to  a  detector  which  may  be  a  photomultiplier  for  detection 
of  wavelengths  in  the  visible  and  near  IR  or  a  photovoltaic  detector 
for  near  and  far  IR  photons.  At  the  same  time,  an  array  of  detectors  is 
irradiated  which  allows  determination  of  misalignment  between  laser  and 
receiver.  This  information  is  sent  to  a  processor  which  issues  com¬ 
mands  to  activators  to  realign  the  two.  The  parameters  e  and  rj  deal  with 
the  losses  through  the  receiver.  These  parameters  are.  functions  of  wave¬ 
length;  thus  they  will  be  computed  for  each  molecule  being  investigated. 
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IV .  Molecules  Investigated 


A.  Sodium 

The  sodium  atom  is  of  interest  in  this  thesis  for  two  reasons. 
First,  it  represents  an  interesting  molecule  for  scientific  reasons. 
There  is  more  natural  sodium  in  the  atmosphere  than  expected.  The  ques¬ 
tion  is,  where  does  it  come  from?  One  theory  is  that  it  is  coming  from 
the  salt  water  in  the  oceans  and  being  drawn  into  the  air  by  thunder¬ 
storm  updrafts  (ref  5:69).  LIDAR  soundings  of  thunderstorms  should  be 
able  to  detect  this  sodium  thus  proving  or  disproving  the  theory. 
Secondly,  it  allows  a  check  of  the  results  of  this  work  with  previous 
LIDAR  sodium  research  (ref  6). 

Transition 

o 

The  (Dj)  (5396A)  line  of  sodium  was  chosen  as  the  wavelength  to 

2  2 

excite  and  detect  the  atoms.  This  electronic  transition  (3  S,  •>  3  Pi) 

2  $ 

is  in  the  visible  (yellow)  region  of  the  spectrum  and  well  within  the 
capability  of  the  dye  laser. 

Effective  Cross-Section 

Assuming  a  sodium  gas  temperature  of  300°K  we  can  determine  the 
doppler-broadened  linewidth  using  equation  (27)  and  the  particular 
values  for  sodium: 

T  =  300°K 

M  =  23  (AMU)  X  t .67  X  lO-27  kg/AMU  =  3.834  X  10"25  kg 

V  -  c_  =  2-9.98  X...IP--CH'/-S.GC  a  5.0843  X  10 14  Hz 
^•0  5896  x  10  cm 
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we  obtain  A\>  =  (2.7341  X  10  IIz) .  We  choose  a  typical  laser  line- 

o 

width  of  1A  (ref  7:350)  which  can  be  converted  into  frequency  band¬ 
width  for  comparison  with  the  sodium  absorption  linewidth  (Av^)  by  the 
formula: 


This  gives: 


Av  =  — rAA 
X1 


AAl  =  8.6242  X  1010  Hz 


(45) 


(46) 


Since  AVjj»Avd  we  may  use  equations  (25)  and  (30)  to  obtain  an  effective 


cross-section  of: 


-17  2 

a  =  1.01  X  10  m 


(47) 


where  the  oscillator  strength  (f)  is  equal  to  (.33)  (ref  26). 


LIDAR  Receiver  Calculations 

Assuming  a  LIDAR  receiver  design  as  in  figure  (4),  we  have  reflect¬ 
ance  losses  from  four  mirrors,  transmission  losses  through  the  collimat- 

o 

ing  lens,  etalon  and  photomultiplier  tube  window.  At  5896A,  a  silver 
coating  appears  to  be  the  best  reflector  for  the  mirrors  (ref  4:187)  with 
a  reflectance  of  (.98).  The  Fabry-Perot  can  be  made  99%  transmissive 
about  a  bandwidth  (AX)  of  1  angstrom.  The  collimating  lens  will  normal¬ 
ly  lose  4%  at  each  surface,  so  it  should  be  about  92%  transmissive. 

The  RCA  7265  photomultiplier  tube  has  a  photocathode  quantum  efficiency 
o 

of  (.09)  at  5896A  ( (Cs)KNa2Sb) .  An  entrance  window  on  the  photomultiplier 
tube  would  have  about  a  92%  transmission  factor.  The  total  throughput 
factor  is  then: 

en  =  (.98)4(.99)2(.92)2(.09)  =  .06885  (48) 

The  diameter  of  the  telescope  primary  mirror  is  initially  chosen  to  be 

one  meter  in  agreement  xvith  several.  LIDAR  design  studies.  This  gives: 
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A  -  .785  m2 

for  the  area  of  the  receiver  aperature. 

Fluorescence  Efficiency 

For  the  sodium  Dj  line,  the  fluorescence  efficiency  (0^)  is  equal 

2 

to  one.  In  other  words,  when  the  sodium  atom  is  excited  up  to  the  3  Pi 

2 

2 

level,  it  can  only  decay  back  to  ground  state  (3  S^)  and  no  other. 
Shuttle  Altitude 

The  shuttle  altitude  is  initially  chosen  to  be  200km. 

Laser  Energy 

The  laser  energy  per  pulse  (I)  is  assumed  to  be  one  joule.  This 
is  an  optimistic  output  energy  for  a  dye  laser  of  about  an  order  of 
magnitude  (ref  7:350). 

Transmission  Coefficients 

Values  of  T  and  T  as  a  function  of  penetration  depth  (y)  into  the 
r  o 

atmosphere  are  generated  by  interpolation  between  values  for  5 145A 
o 

and  6328A,  (ref  1:17,18).  In  this  particular  case,  the  laser  and 
fluorescence  wavelengths  are  the  same,  therefore: 

Tr  =  Tt  =■  T 

Table  II  gives  the  transmission  values  for  the  sodium  Dj  line.  The 
penetration  depth  (X)  is  defined  as  equal  to  zero  at  100KM  altitude  and 
100KM  at  ground  level.  As  an  example  of  how  to  use  the  chart,  if  one 
wanted  to  probe  a  region  70  kilometers  above  the  earth  from  a  shuttle 
platform,  one  would  go  to  Table  II  and  read  off  the  value  for  a  penetra¬ 
tion  depth  into  the  atmosphere  of  30  kilometers  and  thus  T  =  .9991. 

Thus  the  total  transmission  would  be  (.9991)  =  (.9982). 
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Table  II 


Sodium  Dj  Line  Transmission  vs.  Penetration  Depth 


x<M. 

T 

\(KM) 

T 

0 

1.0 

84 

.9779 

30 

.9991 

85 

.974  1 

50 

.9997 

86 

.9449 

55 

.9996 

87 

.9406 

60 

.9994 

88 

.9358 

65 

.9991 

89 

.9307 

70 

.9983 

90 

.9254 

75 

.9961 

91 

.9197 

76 

.9953 

92 

.9135 

77 

.9944 

93 

.9069 

78 

.9933 

9* 

.8971 

79 

.9919 

95 

.8891 

80 

.9902 

96 

.8812 

81 

.9878 

97 

.8571 

82 

.9848 

98 

.8026 

83 

.  .9815 

99 

.6781 

100 

.4323 

Noise  Analysis 

By  utilizing  equation  (41),  the  minimum  number  of  signal  photons 

needed  for  a  given  signal  to  noise  ratio  (R)  can  be  obtained.  For  this 

value,  equation  (9)  can  be  used  to  determine  the  minimum  column  density. 

Before  using  these  equations,  values  for  and,  in  the  sodium 

case,  airglow  (NA),  must  be  calculated.  From  reference  (4:42),  radiance 

o 

values  for  reflected  sunlight  and  moonlight  off  of  the  earth  at  5896A 


are  given  as: 


1.8182  X  1017 


^  photons  ^ 
m7str.&  sec) 


and 

=  2.4242  X  1011  (-j^ons - ) 

*l  m  str.X  sec 


fie*— - 
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The  abbreviation  (str.)  stands  for  seradian.  Sodium  airglow  in  the  atmos¬ 
phere  has  a  daylight  emittance  value  of  20  kilo-Rayleighs^  which  is 
13 

20  X  10  photons  per  square  meter,  and  a  nighttime  emittance  of  100 
12 

Rayleighs  (10  photons  per  square  meter  per  second)  (ref  8).  The 
.  2 

radiance  (photons  emitted  for  m  per  sec  per  steradian)  is  given  by  the 
relationship: 

1  ■  I  <49> 

where: 

M  =  emittance  (20  kilo-Rayleighs  for  sodium  dayglow  and  100  Rayleighs 
for  nightglow) 

The  radiance  (L)  of  the  dayglow  is  then  1,592  X  lo'2  (£~£2H2 - ) 

ci  str.  sec 

and  the  radiance  for  nightglow  is  7.958  X  10 (— ■) .  The  number 

m  str . 

of  photons  detected  per  sample  time  for  each  source  of  noise  is  found 

by  plugging  in  the  appropriate  values  in  equation  (40),  In  the  case  of 

N„  and  N.,  we  have  an  additional  factor  due  to  the  solar  Fraunhofer  line 
E  M 

o 

at  5896A.  This  factor  reduces  and  N,.  because  less  solar  energy  is 
emitted  at  the  Dj  wavelength  due  to  absorption  in  the  sun’s  atmosphere. 
Reflected  moonlight  is  affected  because  that  light  ultimately  comes 
from  the  sun.  Using  equation  (40)  with: 


A  =  .785 

M2 

E  «*  .765 

-4 

>D  "  10 

rad 

AX  =  1A 

n  =  .09 

At  =  6.6  psec 

6  =  3.38 


^ One  Rayleigh  is  defined  as  10 photons  per  tn2  per  second. 


and  a  20%  reduction  due  to  Fraunhofer  absorption  for  N„  and  N  ,  v;e 

E  M 

obtain: 

Nj_,  =  101.88  photons/sample  time 
-4 

N  =  1.36  X  10  photons/sampie  time 
-2 

=  4.46  X  10  photons/sample  time  (day)  and 
7.10  X  10  ^  photons/sample  time  (night) 

The  dark  current  noises  (N^)  can  be  obtained  by  evaluating  the  Richard¬ 
son  equation  (34)  and  substituting  the  result  into  equation  (35).  The 
work  function  ( )  of  the  RCA  7265  photocathode  was  obtained  by  determin¬ 
ing  where  the  response  goes  to  zero  on  the  response  vs.  wavelength 

o 

curve  (ref  10).  This  occurs  at  8200A.  By  using  the  relationship: 

\|>(ev)  =  (50) 

\(A) 

o 

with  X  =  8200A,  one  obtains  1.51  electron  volts  for  \jj.  Assuming  a 

o  2 

photocathode  cooled  to  223  K  and  a  photocathode  area  of  (.785)  cni 

(1  cm.  in  diameter),  we  obtain  1.92  X  10  electrons  per  sample  time 

for  N^.  Therefore,  dark  current  can  be  neglected.  During  daylight 

observations,  it  is  clear  that  Nr,  dominates  all  other  sources  of  noise. 

We  may  therefore  neglect  contributions  of  and  NA  in  equation  (41). 

Using  equation  (41),  a  table  of  N  vs.  R  can  be  generated  for  daylight 

observations  (Table  IV).  For-  night  LIDAR  soundings,  Ng  =  0  and  Nf(  and 

become  the  dominant  noise  sources.  Dark  current  again  can  be  ignored. 

Using  equation  (41)  again,  a  table  of  N  ^  vs.  R  for  night  observations 

can  be  generated  (Table  V).  Choosing  R  =  10,  tables  of  minimum  detectable 

column  density  (K)  vs.  range,  look  angle,  laser  power  and  atmospheric 

penetration  (Tables  VI  thru  IX),  are  presented. 
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Table  IV 


Daylight  N  .  vs .  R  -  Sodium  Case 
J  “  min 


• 

min 

R 

3.88 

1 

7.90 

2 

16.39 

4 

45.80 

10 

109.61 

20 

149.32 

25 

Table  V 

Night  N  .  vs.  R  -  Sodium  Case 
6  min 


11  . 
mm 

R 

.  181 

1 

.568 

2 

2.227 

4 

14.2 

10 

56.8 

20 

88.76 

25 
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Table  VI 


Daylight  Minimum  Detectable  Column  Density  vs.  Laser 
Energy  -  Sodium 


Energy/Pulse 

(Joules) 

I 

,jq!6  photons, 
pulse 

K  . 
min 

( io12m-2) 

n  . 
min 

( 103  cm~3) 

1 

297 

3.6 

3.6 

.8 

238 

4.5 

4.5 

.6 

178 

6.0 

6.0 

.4 

1 19 

9.0 

9.0 

.2 

59 

18.2 

18.2 

.  1 

30 

35.9 

35.9 

.01 

3 

359.0 

359.0 

.002 

.6 

1790.0 

1790.0 

.001 

.3 

3590.0 

3590.0 

At  = 

6.6  ysec 

T2  = 

.98 

AZ  = 

1KM 

A  = 

.785  M2 

i  7 

n  = 

.09 

O  = 

1.01  x  io-' 

Z  =  100KM  e  =  .765 

0=0°  (nadir)  Shuttle  altitude  =  200KM 

R  =  10 

Table  VII 


Nightime  Minimum  Detectable  Column  Density  vs.  Laser 
Energy  -  Sodium 


Energy/Pulse 

(Joules) 

I 

j.  jq  16  photons ^ 
pulse 

Kmin 

(io'2m~2) 

n  . 
mm 

(  103  cm 

1 

297 

1.1 

1.12 

.8 

238 

1.4 

1 .40 

.6 

178 

1.9 

1.87 

.4 

1 19 

2.3 

2.80 

.2 

59 

5.7 

5.65 

.  3 

30 

1 1.3 

1  1 .3 

.01 

3 

i  13.0 

1  13.0 

.002 

.6 

565.0 

565.0 

.001 

.3 

1130.0 

1 130.0 

Other  parameters  are  the  same  as  those  for  Table  VI. 

34 


Table  VIII 


Daylight  Minimum  Detectable  Column  Density 


vs.  Range  and  Look  Angle  -  Sodium  Case 


z 

Z 

V 

K  . 
min 

2 

T  secQ 

0 

n  . 
nun 

(KM) 

(KM) 

( io13m-2) 

(Deg) 

( 10^ cm  3) 

200 

200 

7.7 

.  185 

0 

7.7 

198 

198 

2.2 

.640 

0 

2.2 

196 

196 

1 .8 

.774 

0 

1.8 

194 

194 

1.7 

.810 

0 

1.7 

187 

187 

1.4 

.884 

0 

1.4 

175 

175 

1.1 

.980 

0 

1.1 

150 

.82 

.980 

0 

.82 

125 

125 

.57 

.980 

0 

.57 

100 

.36 

.980 

0 

.36 

75 

75 

.20 

1.0 

0 

.20 

230.94 

200 

13.3 

.  142 

30 

13.3 

228.63 

198 

3.1 

.597 

30 

3.1 

226.32 

196 

2.4 

.744 

30 

2.4 

224.01 

194 

2.3 

.784 

30 

2.3 

215.93 

187 

1.9 

.867 

30 

1.9 

202.07 

175 

1.5 

.977 

30 

1.5 

173.21 

150 

1.1 

.977 

30 

1.1 

144.34 

125 

.76 

.977 

30 

.76 

115.47 

100 

.49 

.977 

30 

.49 

86.60 

75 

.27 

1.0 

30 

.27 

282.84 

200 

31.0 

.092 

45 

31.0 

280.01 

198 

5.2 

.532 

45 

5.2 

277.18 

196 

3.9 

.697 

45 

3.9 

274.36 

194 

3.6 

.742 

45 

3.6 

264.46 

187 

3.0 

.339 

45 

3.0 

247.49 

175 

2.2 

.972 

45 

2.2 

212.13 

150 

1.6 

.972 

45 

1.6 

176.78 

125 

1.1 

.972 

45 

1.1 

14  1.42 

100 

.73 

.972 

45 

.73 

106.07 

75 

.40 

1.0 

45 

.40 

I  =  1  joule/pulse 

All  other  fixed  parameters  are  the  same  as  those  at  the  bottom  of  Table  VI 


Table  IX 


Nightime  Minimum  Detectable  Column  Density 


vs.  Range  and  Look  Angle  -  Sodium  Case 


z 

Z 

V 

K  . 
min 

T2sec0 

Q 

n  . 
min 

(KM) 

(KM) 

( io12m-2) 

(Deg) 

( lO^cm 

200 

200 

23.8 

.  185 

0 

23.3 

198 

198 

6.8 

.640 

0 

6.8 

196 

196 

5.5 

.774 

0 

5.5 

194 

194 

5.1 

.310 

0 

5.1 

187 

187 

4.3 

.884 

0 

4.3 

175 

175 

3.4 

.980 

0 

3.4 

150 

150 

2.5 

.980 

0 

2.5 

125 

125 

1.8 

.980 

0 

1.8 

100 

100 

1.1 

.980 

0 

1.1 

75 

75 

6.2 

1.0 

0 

6.2 

230.94 

200 

41.2 

.  142 

0 

41.2 

228.63 

198 

9.6 

.597 

0 

9.6 

226.32 

196 

7.6 

.744 

0 

7.6 

224.01 

194 

7.0 

.784 

30 

7.0 

215.93 

187 

5.9 

.867 

30 

5.9 

202.07 

175 

4.6 

.977 

30 

4.6 

173.21 

150 

3.4 

.977 

30 

3.4 

144.34 

125 

2.4 

.977 

30 

2.4 

115.47 

100 

1.5 

.977 

30 

1.5 

86.60 

75 

8.3 

1.0 

30 

8.3 

282.84 

200 

95.8 

.092 

45 

95.8 

280.01 

198 

16.2 

.532 

45 

16.2 

277.18 

196 

12.1 

.697 

45 

12. 1 

274.36 

194 

11.2 

.742 

45 

11.2 

264.46 

187 

9.2 

.839 

45 

9.2 

247.49 

175 

6.9 

.972 

45 

6.9 

212.13 

150 

5.  1 

.972 

45 

5.1 

176.78 

125 

3.5 

.972 

45 

3.5 

141.42 

100 

2.3 

.972 

45 

2.3 

106.07 

75 

1.2 

1.0 

45 

1.2 

All  other  parameters,  same  as  Table  VIII. 

The  minimum  column  density  values  given  above  are  compatible  with 

13  -2 

obtained  in  ref  6.  They  detected  a  column  density  ol  6  X  10  M 


similar  LIDAR  equipment. 


results 

using 


CW  vs.  Pulsed  LIDAR  Mode 


An  analysis  was  made  to  determine  for  a  given  scenario,  which  laser 
mode,  CW  (continuous  wave)  or  pulsed  gave  the  larger  signal  to  noise 
(R)  thus  allowing  detection  of  a  smaller  column  density.  The  scenario 
involves  a  one  second  sounding  of  the  atmosphere  around  nadir  from  a 
shuttle  altitude  of  200KM.  A  of  one  second  was  again  chosen  to  allow 
holding  0  constant;  thus  a  potentially  complicated  integration  over  angle 
is  eliminated.  From  the  data  of  Bowman  (ref  6)  it  was  determined  that 
the  bulk  of  natural  sodium  in  the  atmosphere  lies  between  80  -  100  KM 
altitude . 

The  total  number  of  signal  photons  obtained  in  the  CW  case  is  given 
by  equation  (1)  where  I  now  becomes  Iqt^  where: 

I  =  laser  output  power  ( 1  watt) 

=  total  sounding  time  ( 1  second) 

and  the  integration  of  K  is  between  80  and  100  KM.  We  make  the  approx¬ 
imation  that  Z  =  100  KM  (dsitance  from  the  shuttle  to  the  midpoint  of 
the  cloud).  From  ref  6,  the  column  density  is  found  to  be  6  X  10^1 
Substituting  in  the  values  of  K  and  Z  and  the  other  standard  values  used  in 
Table  VI  we  obtain: 

Ncw  =  640  photons/sec. 

For  the  pulsed  mode,  equation  (1)  is  also  used  in  its  normal  form. 

Choosing  I  =  .1  joule  per  pulse,  f  =  10  Hz  and  At  =>  133  psee  to  allow 
sampling  of  a  20  KM  layer  (equation  10),  we  obtain: 

^pulsed  =  640  photons/total  sounding  time 
So,  for  the  same  amount  of  laser  photons  in  both  CW  and  pulsed  cases, 
we  obtain  the  same  signal  return.  The  signal-to-noiae  however  will  not 
be  the  same.  During  dayLight  observations,  the  number  of  noise 


photons  detected  by  the  receiver  is  given  by  equation  (40).  The  para¬ 
meters  are  all  the  same  as  those  used  previously  except  that  At  is  now 
equal  to  (1  second)  for  the  CW  case  and  133  psec  for  the  pulsed 

case.  This  means  that  the  noise  level  for  CW  will  increase  tremendously. 

Using  equation  (40)  we  have  for  CW: 

N£  =  1.57  X  107  photons 

N^j  =  20.6  photons 

Na  =  6757.6  photons  (dayglow) 

Nl  =  9445  photons 

The  last  term  deals  with  the  additional  noise  due  to  laser  beam 
scattering  off  of  the  earth.  From  (ref  9)  the  irradience  at  the  aperture 
of  the  receiver  due  to  the  emittance  of  a  disk  of  luminance  (L)  of  radius 

(a)  at  a  distance  (Z)  away  where  Z»  a  is  given  by: 


where: 

E  =  irradiance  )  at  the  receiver  plane 

IT 

M  =  emittance  of  the  disk  (jjlgj—— .) 

H 


Recall  that  for  a  Lambertian  source  ttL  =  M 

The  size  of  the  disk  is  determined  by  the  FOV  of  the  detector  (equation 
43).  With  Z  =  200  KM  and  Gq  =  10  4  radians,  we  obtain  a  10  meter  spot 
radius.  The  emittance  (M)  of  the  disk  is  given  by: 

M  =  I  T. r  =  1.6257  X  1014  photons 
o  t 
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where: 


18  w 
I  =  1  watt  *  2.97  X  10  photons/sec  @  5S96A 

T  =  laser  transmission  loss  from  space  to  ground  (.43) 

r  =  earth  albedo  or  reflectance  coefficient  (.04) 

2 

ira  =  area  of  disk  being  illuminated. 

The  number  of.  photons  entering  the  receiver  aperature  per  second  is 
given  by: 

N  =  E  A  =  5.4S76  X  10*^  photons/sec 

where: 

Tr  =  photon  transmission  from  the  earth  back  to  the  receiver  (.43). 
Finally  the  number  actually  detected  in  the  one  second  sounding  time  is 
given  by: 

*  Ne  n  At  =  9445  photons 

It  is  clear  that  the  signal  to  noise  during  daylight  observations  is 
given  by: 

R  =_N_  (52) 

<V* 

due  to  the  overwhelming  magnitude  of  earth  reflected  sunlight.  This 
gives  a  signal  to  noise  of  (.163).  For  the  pulsed  case,  At  =  133  ysec, 
and  thus  the  noise  values  will  be  considerably  reduced: 

Ne  =  -2037.6 

Nm  =  .027 

=  .899  (dayglow) 


Again,  for  daylight  observations  the  signal  to  noise  is  given  by 
equation  (52): 


pulsed 


640 

(2037.6) * 
39 


14. 18 


So  in  daylight  work,  the  pulsed  mode  is  the  only  mode  that  will  work. 
For  night  observations  (N  =0)  we  have: 


640 


cw 


(N+Nl)'2  (640  +  94 4  5)2 


=  6.37 


The  pulsed  case  yields: 


R  ,  .  — V  =  25.3 

puised  (N)^ 


The  pulsed  inode  is  still  the  better  mode  to  use.  It  is  clear  that  for 
nanosecond  radiative  lifetimes  that  the  pulsed  mode  is  far  superior  to 
the  CW  mode  because  all  signal  photons  can  be  collected  quickly  thus  cut¬ 
ting  down  on  noise.  For  long  radiative  lifetimes,  one  may  be  forced  to 
sample  longer  to  get  as  many  signal  photons  as  possible,  but  this  can 
probably  only  be  done  in^ a  noise  free  environment. 


Beam  Divergence 

-4 

A  beam  divergence  of  10  rad  was  assumed  for  the  laser  and  FOV  of 
the  receiver.  The  question  now  is:  "Does  the  laser  need  expanding 
optics  to  reduce  diffraction  and  thus  beam  divergence?"  Using 
equations  (43)  and  (44)  we  obtain  a  required  exit  aperature  for  the  beam 
of  7.19  millimeters.  Such  a  small  size  suggests  that  no  expanding  optics 
are  necessary. 


Sodium  Conclusions 

It  appears  that  very  small  quantities  of  sodium  (average  number 
4  3 

density  of  10  particles  per  cm  )  can  be  detected  with  good  signal  to 
noise  over  large  ranges  and  look  angles.  The  main  reason  for  sodium's 
easy  detection  is  its  large  cross-scct ion  (10  ^).  This  situation 

unfortunately  is  the  exception  rather  than  the  rule.  Therefore  minimum 
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.  .*.•>  *- 


detectable  column  densities  for  other  particles  will  be  higher,  and  in 
some  cases  much  higher  than  sodium.  It  is  also  clear  that  pulsed 
sounding  of  the  atmosphere  is  superior  to  CW  due  chiefly  to  the  large 
reduction  in  noise  in  the  pulsed  mode. 
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B.  Hydrogen  Chloride  (HC1) 

Hydrogen  chloride  is  of  interest  because  it  is  a  possible  by-product 
of  rocket  fuel  combustion.  The  detection  of  certain  molecules  in  a 
rocket  plume  makes  it  possible  to  determine  what  kinds  of  fuels  and 
oxidizers  are  being  used  in  the  rocket  and  this  may  allow  one  to 
ultimately  predict  the  performance  characteristics  of  the  vehicle. 

Transition 

For  a  LIDAR  sounding,  excitation  of  an  electronic  transition  would 
be  the  ■  most  desirable  way  to  detect  HC1  because  of  the  large  associated 
cross-sections.  Unfortunately,  the  excitation  of  an  electronic  transi¬ 
tion  in  HC1  cannot  be  accomplished,  because  the  transition  from  the 
electronic  ground  state  to  the  first  excited  state  lies  far  out  in  the 
ultraviolet  region  of  the  spectrum.  This  is  beyond  the  capabilities  of 
a  tunable  laser.  Thus  one  must  be  content  in  attempting  to  excite 
infrared  vibrational-rotational  levels.  Excitation  of  the  fundamental 
( v"=0->  v,=  l)  and  first  overtone  (v"=0->v’=2)  vibrational  bands  cannot 
be  done  with  a  dye  laser  since  these  transitions  lie  too  far  into  the 
IR,  beyond  1.2  pm  which  is  the  long  wavelength  cutoff  of  the  present 
dye  laser.  The  second  overtone  vibrational  transitions  (v"*0-»-  v  '=3) 
does  allow  use  of  the  dye  laser  because  this  transition  occurs  in 
the  region  of  1.2  pm.  There  are  other  lasers  such  as  the  optical 
parametric  oscillators  ( 0P0)  (ref  1  1:257)  that  allow  tunability  out  to 
3.5  pm,,  however  they  are  a  factor  of  one  thousand  smaller  in  pulsed 
output  energy  and  are  not  as  well  developed  as  the  dye  laser. 

Rotational  transitions  for  the  ground  and  excited  states  must  be 
selected  next.  It  is  desirable  to  excite  the  rotational  level  that  has 


the  highest  population  at  a  temperature  (T) .  This  is  given  by,  (ref  12:41) 

j  -  (— k.T  s'*  _  _L  (53) 

MAX  4hcU;  2 


where : 

k  =  Boltzman's  constant  =  1.38  X  10  erg/°K 

-27 

h  =  Planck's  constant  =  6.626  X  10  erg/sec 
c  =  speed  of  light  =  2.998  X  10 cm/sec 
B  =  rotational  constantf  =  10.5  cm  * 

Using  a  typical  rocket  plume  temperature  of  1500°K  (ref  13:352)  and  a  B 
coefficient  of  10.5  cm  *  for  the  HC1  (ref  12:63),  we  have  for  equation  (53): 


JMAX 


-  7 


Therefore  the  transition  to  be  excited  is: 

(v"=0,j"=7  ->  v'=3,J'=8)tt. -  This  is  the  R(7)  line  of  the  second  over¬ 
tone  transition  at  oj  =  8447.6582  cm  *  or  X  =  1.1838  pm.  It  is  desirable 
to  detect  a  slightly  different  wavelength . so  that  Mie  scattering  by 
particulates  in  the  atmosphere  and  plume  can  be  avoided.  The  transitions 
(vf=3->  vr,=2)  would  ideally  be  the  most  desirable  vibrational  transition 
to  detect.  It  is  in  the  far  IR,  but  this  presents  no  problem  since  many 
detectors  work  quite  well  in  the  IR.  The  problem  lies  in  the  lack  of 
data  on  v'=3-+v"=2  (hot  band)  transitions  especially  data  that  would 
allow  one  to  calculate  an  effective  cross-section.  Therefore,  even 
though  Av=l  or  2  transitions  are  more  likely  than  Av=3,  the  latter  is 


2 

+The  rotational  constant  is  defined  as:  h/8n  Ic  where  I  is  the  moment  of 
inertia  of  the  molecule. 

ft  Note  the  accepted  notation  of  double  prime  refering  to  the  lower 
state  and  single  prime  for  the  upper  state. 
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chosen  simply  because  data  is  available  on  this  transition.  The  rotation 
al  transition  (J'=8->-  J"=9)  at  X  =  1.2202  ym  must  be  used  because  of  the 
selection  rule  J=_f  1  and  because  of  the  desire  to  transmit  and  receive 
at  different  wavelengths.  All  of  the  HC1  molecules  do  not  lie  in  the 
same  excited  state  at  a  finite  temperature.  On  the  contrary,  at  rocket 
plume  temperatures  there  are  many  occupied  states.  The  higher  the  temp¬ 
erature  the  more  states  are  occupied.  It  becomes  important  to  know, 
therefore,  how  many  molecules  there  are  in  the  state  that  is  to  be 
excited  relative  to  the  total  number  of  IIC1  molecules  being  irradiated. 
This  requires  the  use  of  the  Boltzman  distribution  and  an  evaluation 
of  the  partition  functions  for  vibration  and  rotation.  The  Boltzman 
distribution  may  not  be  valid  in  a  rocket  plume  since  the  exhaust  gases 
are  not  in  equilibruim.  However,  the  only  alternative  is  detailed 
accounting  of  all  chemical  reactions  and  many  are  not  well  understood. 
According  to  the  Boltzman  distribution,  the  ratio  of  the  population 
of  the  ith  state  to  that  of  the  ground  state  (nQ)  is  given  as: 

n.  g. 

—  »  _i  exp  (-AE./KT)  (54) 

n  g  i 

o^o 

In  the  case  of  rotational  levels: 

Aeo  =  J(J+l)BhC  (55) 

where: 

J  =  rotational  quantum  number 

=  degeneracy  of  the  ith  level  =  (2i+l) 

B  =*  rotational  constant  (for  HC1)  =  10.59  cm  *) 

The  ratio  of  interest  is  nj/nTQT  where  represents  the  total  number 

of  molecules.  The  total  number  of  molecules  can  be  expressed  in  terms 
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of  by  using  the  Boltzinan  distribution: 

n„,„„  =  n  +n,+n„+. .  .+n.  +  « . . 
10T  o  1  2  i 


®  i 

ni  =  no  7"  exp  (  -aeAD 


n  g . 

nTOT  “  no  +  ~"V  ~~  exp  + - + 


n  g- 
o  &i 


exp  (-AE./KT) 


no®N 

+ _ +  exp  (  Af^/KT) 

so 

where  the  Nth  level  is  the  one  of  highest  energy. 
One  can  simplify  the  above  equation  to  the  form: 


TOT  '  g 
eo 

where  Z  is  known  as  the  partition  function: 


Z  =  E  gi  exp  (-AE./KT) 
i=o  1 

Returning  to  equations  (46)  and  (43),  we  have: 

ntot  sniF  exp  (-^i/KT)  f- 

eo 

which  yields: 

n . 

=  g.  exp  (-AE./KT)  (57) 

TOT  - ^ - - - 

The  ratio  of  the  number  of  molecules  in  the  J'th  rotational  level  relative 
to  the  total  number  of  molecules  in  the  v=0  state  is  obtained  using 
equations  (55)  and  (57) 


*  exp  (-J ' (J *+  t )RhC/KT)  (2J'+1)  (58) 


where : 


(59) 


N 

ZR  =  I  (2J+  1 )  exp  (-J(J+t)Bch/KT) 

J=0 

Using  the  same  reasoning  for  vibrational  transitions,  we  have  initially: 

AE^  =  hcvue  (60) 

where: 

V  =  vibrational  quantum  number 

to  =  vibrational  constantf  =  2.990.9463  cm  '  (for  HC1) 
e 

In  the  same  form  as  equation  (57),  we  have  for  vibration:  • 


nv 1  =  exp  (“hcv  toe/KT)  (61) 

nT0T  Z 

V 

Here  nT0T  represents  the  total  number  of  molecules  in  all  vibrational 
levels , 

The  vibrational  partition  function  is  given  by: 

N 

Z  =  Z  exp  (-hcv  u)  /KT)  (62) 

v=0 

Vibrational  levels  are  not  degenerate;  therefore  the  g^  terms  do  not 
appear.  It  is  clear  that  by  multiplying  equations  (58)  and  (61)  to¬ 
gether,  the  ratio  of  the  number  of  molecules  in  the  (v',J')  state  to  the 
total  number  of  HCl  molecules  is  given  by: 


V,J’ 


=  (2J'+1)  exp  (-J 1 ( J 1 + 1 )Bhc/KT)  exp  (-hcu'Up/KT) 


tue  is  known  as  the  hypothetical  equilibrium  frequency  of  the  anharmonic 
oscillator.  It  is  the  frequency  for  infinitely  small  vibrations  about 
the  equilibrium  point  of  the  molecule. 
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The  evaluation  of  the  partition  functions  Z  and  Z  can  at  times  become 

K  V 

very  tedious  unless  some  approximations  are  made.  This  can  be  done 
for  the  vibrational  case  but  not  for  rotation.  For  vibration  at  1500°K, 
it  can  be  shown  that  most  of  the  molecules  are  in  the  v=0  state.  Using 
the  Boltzman  distribution  and  equation  (54)  we  have: 


and 


=  exP(— |2*L)  =  .056684 
v=0  iv 


n 

v= 


2 


v=0 


exp  (• 


-2hcioe 

KT 


)  - 


,0032 


From  these  ratios,  it  is  clear  that: 


V-U  1  10i 

Using  equation  (53),  we  have 


nv=0  ~  exp  (-(0)hcuja/KT) 

1,101  E  exp  <-hcv5  /KT) 

v=0  e 


~  _ 1 _  =  .946 

nTQT  1+. 056684 

In  the  case  of  rotation  at  1500°K,  many  levels  are  substantially  pop¬ 
ulated  and  no  convenient  approximation  can  be  made.  Evaluation  of  the 
first  25  terms  of  the  partition  function  was  considered  sufficient 
accuracy  for  this  study  because  higher  energy  levels  are  thinly  pop¬ 
ulated  and  contribute  a  negligible  amount,  to  the  summation. 

It  is  found  that  for  HC1: 

25 

ZD  ~  l  (  )  :  60.576 


I 
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and : 


nTOT 

V=0 


(2(7)+ I )  t:xp  (-2(7+  OHpC./KT) 
60.5/6 


.  141 


Finally,  the  ratio  of  n  to  n  is  given  as: 

V  —  U ) J - <  Awl 


-v"  ° ~7  =  (.946)  (.141)  =  .  133 

nT0T 

or: 


nv=0,J=7~n 


.  133  n 


TOT 


Effective  Cross  Section 

The  absorption  cross-section  is  the  next  item  that  must  be 
calculated.  Unfortunately,  the  Einstein  coefficients  were  not  avail¬ 
able  for  HC1  to  allow  a  quick  calculation  of  the  effective  cross-section. 
However,  data  on  dipole  transition  matrix  elements  for  the  fundamental 
and  first  and  second  overtone  transitions  were  given  for  11C1  and  these 
could  be  used  to  calculate  an  A'  coefficient  and  then,  finally,  an 
absorption  cross-section.  The  frequency  of  a  spontaneous  transition  from 
an  upper  level  (v',J')  to  a  lower  level  (v",J")  is  given  by  the 
expression:  (ref  14:121). 


v",J" 


a/ 3 

64tt 

=  3h (2J ' V 1 ) 


[”l— 8BM  nf 

ffinr- 

~i e  1  e 


(64) 


where : 

io  =*  band  origin  in  cm  * 
o 

h  =  Planck's  constant  (erg-sec) 
m  =  -J"  (P-branch  transitions) 


or: 

m  =  J"+l  (K-branch  transitions) 
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|  Rv  [  =  rotat  ionless  dipole  moment  transition  matrix  element  (esu-cm) 
v" 

M  =  Dipole  Moment  coefficients 

o’  1  1 

rg  =  equilibrium  distance  between  atoms  of  the  molecule 

The  last  factor  in  equation  (64)  is  sometimes  known  as  the  Herrnan- 

Wallis  correction  factor  for  rotation. 

Before  evaluating  A',  the  dipole  transition  matrix  element  had  to 

be  calculated.  The  dipole  transition  matrix  element  is  defined  as: 

00 

Rv  =  K)*„  fr)  M(r)  \5  <  f r )  dr  (65) 

"  JYv  » J'>  Vv  >J* 

V  o 

Numerical  solution  of  the  Schrodinger  wave  equation  using  a  quintic 
anharmonic  oscillator  potential  and  third  order  perturbation  theory 
was  used  (ref  15)  to  obtain  expressions  for  the  wavefunctions  m  T 

YV,J . 


A  cubic  dipole  moment  function  of  the  form: 

M(r)=Mn+M.(r-r  )+M„(r-r  )2+M,(r-r  )3 
0  1  e  i  e  3  e 

was  used.  The  evaluation  of  the  integral  in  equation  (65)  gives  finally 
for  the  second  overtone  rotationless  transition  matrix  element: 


Y  = 


0. 


(3) 

8 


3/ 2,..  '  3  2. 

f  (Mlre(a2+  4  31  } 


4M0r 
2  e 


4M3re3) 


a  j  =  -2.364 
a2  =  3.662 
a3  =  -4.710 


.00708 


o 

.914  Debye/A 
°2 

.025  Debye/A 

°3 

-1.54  Debye/AJ 


_  j  g 

( 1  Debye  =  10  esu-cm 
This  yields: 

Rq  =  -2.  1559  X  IQ-21  esu-cm 
Using  Equation  (64),  the  Einstein  coefficient  becomes: 

A,(v'^,J'«.8+V"*fl,J".7)  =  .31396  sec"1 

The  effective  absorption  cross-section  is  given  again  by  equations  (26) 
(27),  and  (28).  The  line  shape  profiles  a(^)  and  L(v)  are  assumed  to 
be  reasonably  well  represented  by  a  Doppler-broadened  profile.  The 
width  of  the  HC1  absorption  line  is  given  by  equation  (19)  with: 

M  =  (35AMU)  (1.66Xl0~27)kg/AxMU  =  5 .8  IX  10_26kg 

T  =  1500°K 

v  =  C/X  =  C/ 1.1838  pm  =  2.5325X1014HZ 
o 

we  have: 

AV  HCl  =  1.  1873XI09HZ. 


A  typical  dye  laser  linewidth  of  one  angstrom  is  again  assumed,  or  Av^ 
2 1 .393X109HZ.  We  can  reasonably  say  th3t  Av  ^  »  Av"^  and  make  the 


same  approximation  for  °e££  as  was  done  in  the  sodium  case  by  using 
equation  (30). 

This  yields: 

0  ..  =  8.7093X10"21  cm2  =  8.7X10~25M2 
ef  f 

Note  that  this  cross-section  i.s  a  factor  of  about  ten  million  smaller 
than  that  of  sodium. 
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The  lifetime  of  the  fluorescence  transition  ( v '=3 , J ' -8 +  v"=0, J"-9) , 
which  is  the  P(9)  line  of  t!ie  second  overtone  is  given  by  equation  (64) 
with  m--9: 

A'  =  .55637  sec'1 

f luor . 

This  means  that  the  lifetime  of  the  state  is: 

t  =  tt  =  1  *797  sec 
A 

This  is  a  very  long  lifetime,  so  equation  (14)  should  be  used  to 
determine  the  total  number  of  signal  photons  collected  in  a  sounding 
period  (i^). 

Using  the  parameters: 

1  Q 

I  =  1  joule/pulse  =  5.9617X10  photons  @  1.1838  pm 
At  =  6.6  sec 
N  =  10  pulses 
A’  =  .55637  sec-1 
a  =  8.7093X10-25  M+2 
Z  =  300  KM 
A  =  .785  M2 
ri  =  .  1 
e  =  .765 
0F  -  .05 
Az  =  10  m 

T„  =  T.  =  1  (assuming  rocket  at  altitudes  above  70  KM) 

K  t 

we  have : 

=  (1.3766X10-20)  (  —  )  K  (  E  E  exp  (-JA*/10) 

TOT  n 

Det.  i=0  J=0 
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The  quantities  (n -/nT)  and  A*  can  be  evaluated  by  first  determining  the 

mean  time  between  collisions  (-j  )  by  using  equations  (16)  and  (17). 

From  a  rocket  plume  computer  model  run  by  the  foreign  technology  division 

(FTD),  the  total  number  of  colliding  molecules  (nc>  is  about  (5.8S) 

1 7  -3 

times  larger  than  the  total  number  of  HCl  molecules  (~)0  cm  ).  The 
average  velocity  of  the  HCl  molecules  is  given  by  equation  (16)  with: 
m  =  5 .8 lX10~26kg 
T  =  I500°K 

This  yields  a  velocity  of  1.034X10"*  cm/sec. 

The  approach  now  is  to  choose  a  column  density  (K)  and  solve  for 

Det . 

This  is  simplier  than  choosing  a  and  arriving  at  a  column  density 

Det 

since  much  more  algebra  is  involved.  When  a  value  for  column  density 
is  specified,  a  value  for  HCl  particle  density  (n)  and  total  collision 
particle  density  (n^)  is  also  specified  by  the  following  relationships: 

K  =  n'AZ  =  (.133)  nAZ 

and 


n  =  5.88n 
c 

If  the  plume  is  modeled  to  be  10  meters  thick,  then  &Z  =I0  meters.  With 

the  above  values,  a  table  of  column  density  vs.  is  constructed. 

Det 
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Table  X 


Column  Density  vs.  N 


K 

n 

N 

TOT 

o 

Det 

(m-2) 

(cm  ) 

(counts) 

1.33  X 

io’6 

io10 

1.00  X  10"U 

1.33  X 

,o20 

V4 

5.54  X  10“5 

1.33  X 

io23 

io17 

1.68  X  10-5 

1.33  X 

1025 

1019 

1.68  X  10-5 

From  the  results  in  Table  X,  it  is  clear  that  detection  of  HCl  in  the 
scenario  envisioned  is  impossible  for  reasonable  column  densities.  The 
main  reason  for  the  problem  is  the  small  Einstein  transition  coefficient 
and  the  resultant  small  absorption  cross-section  and  long  radiative  life¬ 
time.  Small  column  densities  result  in  less  quenching  but  the  slow 
radiative  lifetime  of  the  state  results  in  a  small  number  of  radiative 
transitions  occuring  in  the  total  sounding  time  period  (t^).  Large 
column  densities  do  not  help  the  signal  because  collisional  deexcitation 
deexcites  the  molecule  before  it  has  a  chance  to  radiate.  This  result 
manifests  itself  in  the  small  value  of  the  ratio  (nD/n„)  in  the  modified 
LIDAR  equation.  One  could  hope  to  obtain  more  signal  counts  by 
continuously  sampling  the  return.  In  other  words  At -*-f  ' .  Of  course 

in  this  case,  one  loses  spatial  resolution  capability  (AZ)  and  one  is 
also  allowing  much  more  noise  into  the  system.  This  approach  could  only 
be  done  in  a  scenario  almost  totally  free  of  background  noise.  However 
this  still  would  not  help  as  the  signal  counts  would  never  exceed  IO-4 
in  the  one  second  sounding  period  which  can  be  verified  by  substituting 
(.01)  in  for  At  and  solving  the  modified  LIDAR  equation  with  the  same 
parameters  as  before. 
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Conse q nonces  of  a  Laser  With  a  broader  Tunable  Range 

If  a  tunable  laser  were  developed  that  could  reach  into  the  ultra¬ 
violet,  then  it  would  be  possible  to  take  advantage  of  the  laser 
absorption  cross-sections  for  electronic  transitions  and  detection  of 
HCl  through  resonance  fluorescence  would  probably  be  feasible.  If  a 
tunable  laser  is  developed  that  could  look  farther  into  the  IR,  one 
could  take  advantage  of  the  larger  absorption  cross-sections  and 
shorter  radiative  lifetimes  of  the  fundamental  band  of  HCl.  If  one 
excited  the  1-0  band,  the  absorption  cross-section  would  be  increased 


by  a  factor  by: 


OJ0 

°30 


*  9  2  3,2 

A10X210  _  |Ro  1  co10X 

A30X230  "  |R2|2o4  X2q 


378 


where: 

AI0  =  48.3  A30  :  19  sec 


^ 10  ~  3030 


-1 

cm 


Xjq-  3-3  ym 


Also,  for  a  laser  output  of  1  joule,  one  would  obtain  a  factor  of  2.8 
more  photons  at  3.3  ym  than  at  1.18  ym.  Also,  the  quantum  efficiency 
(0  )  would  increase  by  an  order  of  magnitude  since  a  v'=l  level  can  only 

f 

decay  to  the  v"=0  level.  Substituting  these  values  into  the  modified 
LIDAR  equation  with  the  parameters  used  above  gives  about  60  signal 
counts  for  ten  sample  times  6.6  ysec  in  duration.  This  is  indeed  a  good 
signal  and  makes  detection  of  HCl  feasible  if  a  laser  can  be  found  to 

t 

excite  the  (1-0)  band. 


Conclusions 

Detection  of  HCl  through  LIDAR  fluorescence  techniques  appears 
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unlikely  unless  lasers  of  reasonable  output  energies  can  be  built 
to  excite  electronic  transitions  in  the  UV  or  IR  -  vibrational- 
rotational  lines  in  the  fundamental  band.  Present  dye  lasers  can  only 
excite  overtone  bands  of  HC1  and  these  have  absorptions  cross-sections 
too  small  for  a  reasonable  L1DAR  return. 


C.  Hydrogen  Flouri.de  (HF) 


Choice  of  Transit  ion 

Another  molecule  which  may  appear  in  a  rocket  plume  is  HF.  The 
first  task  again  is  to  choose  a  transition  for  excitation  and  fluorescence. 
Unfortunately,  like  HC1,  all  electronic  transitions  in  HF  are  in  the  far- 
ultraviolet  and  not  accessible  to  the  dye  laser.  One  is  forced  again 
into  excitation  of  vibrational-rotational  transitions  and  must  live  with 
the  resultant  smaller  absorption  cross-sections.  Given  a  temperture 
of  1500°K,  equation  (53)  can  be  used  to  determine  the  most  populated 
rotational  level.  In  the  case  of  HF,  B  =  20  cm  *  (ref  12:45)  and  sub¬ 
stituting  this  into  equation  (53)  yields: 


The  number  density  of  nJ5  relative  to  n^^  can  be  found  by  utilizing 

v=0 

equation  (58).  This  gives 


J=5 


(.1)8)  n, 


TOT 

V=0 


The 


ratio  of  n  ,  to  n  _  is  given  by: 
V=  1  v=0  6 


V= I  ~  ,-hCioe^ 

—  -  exp  (  jjy-> 

V=0 


-1 


using  u  =  4138.5  cm  for  HC1  we  have: 
e 


'  % 
1 

1v=0 


.019 


It  is  clear  therefore  that  n^_g  ~( . 98)  n^^  and  the  ground  state  of 
interest  is  therefore  v"=0,  J"=5,  The  ratio  of  ground  state  molecules 
to  the  total  number  of  molecules  is  then,  by  equation  (63),  equal  to: 
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v>"=0  I"  =  5  n' 

-iLJid — d.  „  IL_  „  (.MS)  (.98)  =  (.116) 

nTOT  nTOT 

A  transition  is  needed  such  that  the  exciting  wavelength  will  be  within 
the  capability  of  the  dye  laser.  Excitation  from  the  v"~0  to  \>'=1  or  2 
is  too  far  out  in  the  IR  for  the  laser.  Again,  as  in  the  case  of  HC1, 
the  second  overtone  vibrational  band  must  be  considered  if  one  uses  a  dye 
laser.  Recalling  the  selection  rule:  Aj=+.  1,  the  excitation  transition 
is  chosen  to  be  (v"=0,J"=5 ->■  v' =3 1  =6) .  This  is  the  R(5)  transition. 

Spectral  line  data  on  the  second  overtone  band  was  not  available, 
however,  spectral  lines  could  be  calculated  by  using  the  Dunham  co¬ 
efficients.  For  the  laser  line  we  have  (ref  12:74): 

||  -  Ae  -  j  cm"1  =  (2B(J"+ 1)-4D  (J"+l)3  +  3wg  (l-4Xe)) 

where : 

Y  =  vibration  correction  factor  for  the  anharmonic  oscillator 
e 

j  m  (2(20)(6)-(4)(.00187)(6)3+3(4138.5)(l~4(  .02 18) ))  cm" 1 

|  =  11571.253  cm-1 
Finally: 

o 

Alaser  =  8642 . 1  A 

From  (ref  14:128),  the  (v'  =3  -*■  v"=2  or  (v'=3  -*■  V"=  1 ) 
fluorescence  are  much  more  likely  to  occur  than  the  second  overtone 
transition  since  these  transitions  have  much  higher  Einstein  A’  co¬ 
efficients.  Unfortunately,  the  v'=3-*-v"=2  band  which  has  the  highest 
A*  coefficient  of  the  three  is  not  a  desirable  wavelength  regime  to 
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monitor  because,  like  the  fundamental  band,  it  lies  in  the  middle  of 
the  2.7  ym  water  emission  band  which  could  introduce  unacceptable 
absorption  and  noise  into  the  LIDAR  experiment.  Thus,  the\7=3,  \>''=1 
transition  is  chosen. 

Utilizing  the  selection  rule  AJ=  _+  1  one  may  observe  either  the 
(v'=3,J' =6 ->  ,J"=7)  or  the  (v'=3, J '=6  -*■  v"=  1  ,J"=5)  which  are  the 

P(7)  and  R(5)  transitions  respectively.  The  P(7)  transition  is  arbitrari¬ 
ly  chosen.  The  wavelength  of  this  transition  is  given  by: 

X(cm)  =  (Aev  -Aej)  *  (66) 

where : 

Aev  «*  Energy  difference  between  the  vibrational  levels. 

Ae  *  (v1  +  4-)t0  -Xo(v'  +  t)2w  ”(v"  +  4-)u 

V  2  e  Ae  2  e  2e 

/  H  .  K2 

+Xe(v  +  ~e 

Ae  ,  .  „  .  -  ((3.5)  (4138.5)  -  (.0218)  (3.5)2  (4138.5 

V  C  9  V  =  * 

-  (1.5)  (4138.5)  +  (.0218)  (1.5)2  (4138.5)) 

Ae^»  =3>  vn  =  j  =  7374.807  cm  1 

Calculating  the  rotational  energy  gaps  and  making  use  of  the  fact  that 
the  rotational  correction  constant  is: 

D  . 


we  have: 


A Cj  -  OU"(j"  +i)-DJ"2(J"+l)2-BJ'(j'+i)tyj,2(j’+))2) 


Aej i_.6 ^ j"=7=(^2°) (56)=( .00lS7)(49)(64)-(20)(42)+(.001S7(36)(49)) 


Aej  =  277.43436  cm 

Returning  to  equation  (66)  we  have  for  the  fluorescence  wavelength: 

XF(cm)  =  1.4O9XI0-4  cm  =  1.4090  ym 

This  transition  is  in  the  I.R  and  an  InSb  1R  detector  is  better  to  use 
than  a  photomultiplier  tube  since  the  quantum  efficiency  at  the  fluor¬ 
escence  wavelength  is  approximately  (.27). 


Absorption  Cross-Section 

The  rotationless  Einstein  A’  coefficient  for  the  laser  wavelength 
was  obtained  from  (ref  14:128).  The  coefficient  does  not  include  the 
Herman-Wallis  correction  factor  for  rotation,  however  this  correction 
is  usually  very  small  and  can  be  neglected.  The  value  is  1.223  sec 
Using  equation  (28)  to  calculate  Avu„  at  T  =  1500°K  gives: 

nr 

A\)  =  2.1516X109HZ 

tit 

With  a  laser  linewidth  of  AX=100pm  or  Av  =  ,40141X109HZ  it  is  clear 
that  Av^  and  AvHF  are  close  to  being  equal.  Therefore,  the  complete 
expression  for  the  effective  cross-section,  equation  (21),  should  be 
used.  This  yields: 

-?2  2 

a  ,,  =  1.0139  X  10  ^  M 
err 

The  Einstein  A’  coefficient  for  the  3-1  vibrational  transition  is  equal 
to  65.44  sec  *  (ref  14:128). 


if’ 
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One  can  now  utilize  equation  (14)  to  determine  K  vs,  N, 


Choosing  the 


TOT  ‘ 
Det 


following  parameters: 

I  =  4.3521X1018  Ph~-?— -  @  8642. 1A 
pulse 

°  ,,  =  1.0139X10_22m2 
elf 

A'31  =  65.44  sec"1 

At  =  6.6  v<sec 
A  =  .785  m2 
Z  =  3X105m 
n  =  .27 


£  =  .765 


n '  =  . 11 8n 

N  =10  pulses 

-15  2 

O  =  10  cm 
c 

AZ  =  10m 


n£  *=  2ii  (from  FTD  plume  model) 
in  equation  (14),  with  0^  determined  by: 


A' 
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A.' 3 1+A  30+A  32 


=  .069t 


we  obtain: 

Ntot  =  (1.3095X10-17)  (Us.)  k  O-e  A*At)  X  E  exp  (-JA*/!0) 

Det  nT  i=0  J=0 

Again,  as  in  the  case  of  HC1,  one  assumes  a  column  density  (K)  and 

solves  for  which  is  the  number  of  signal  photons  or  counts  registered. 

Det 


tThe  factor  of  (^)  in  the  expression  for  ££,  arises  because  of  the  two 
equally  probable  rotational  transitions  that  can  occur  for  each  vibrational 
transision  . 
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Table  XI 

HF  Column  density  00  vs.  Signal  Counts  (N 

Det 


K 

n 

ntot 

(m“2) 

(cm  ^) 

Det 

(counts) 

1.18  X  1016 

io10 

6.68  X  10"' 

1. IS  x  io!8 

io12 

5.36  X  10 

1. 18  X  102° 

io1A 

6.10 

72 

1.18  X  10 

10 16 

36.90 

23 

8.73  X  10^ 

7.4  X  1017 

36.90 

1 7  -3 

The  HF  number  density  of  7.4  X  10  cm  that  appears  in  Table  XI  comes 
from  the  FTL)  rocket  plume  model  mentioned  briefly  in  the  HC1  section. 

From  the  above  table,  it  is  clear  that  for  reasonable  HF  number  densities 
(10  cm  ),  LIDAR  detection  of  HF  seems  feasible.  The  next  question 

is:  Is  background  noise  at  the  fluorescence  wavelength  low  enough  to 
achieve  a  signal  to  noise  greater  than  one? 


HF  Noise  Analysis 

From  the  ESA  report  (ref  4:42)  the  radiance  values  for  terrestially 
reflected  sunlight  and  moonlight  at  the  fluorescence  wavelength  (1.4  ym) 
are  given  as : 

Lg  -  1.3636  X  I017  - 5 

m  str.  sec  A 


and : 

=  2.4242  X  1011  F 

m  str.  sec  A 
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Using  equation  (40)  with: 

=  3.33  X  10  ^  rad  (yields  a  spot  size  of  10  meters  at  300km) 
o 

AX  =  1A 

we  can  determine  the  number  of  background  photons  detected  to  be: 

N„  =  127.3  photons  detected/sample  time 

=  2.32  X  10  ^  photons  detected/sample  time 

Using  these  values  in  equation  (41)  without  the  factors  of  £  since  a 
photomultiplier  tube  cannot  be  used  here  at  1.4  ym,  we  can  determine  the 


minimum  number  of  signal  counts  (K^^)  per  sample  time  to  obtain  a  certain 
signal  to  noise  (R).  For  daylight  observations,  dominates  all  noise 
terms  and  at  night  we  have  only  N^.  is  again  considered  to  be 

negligible.  With  the  above  values,  tables  of  N  vs.  R  are  given  for 


day  and  night  soundings. 


Table  XII 


Day  N  .  vs.  R  at  1.4  p 
min  ^ 


(Counts/sample  time)  Signal/Noise) 


11.16 


15.09 


Table  XIII 


Night  N  .  vs.  R  at  1.4  ym 
“  min 


N  . 
mm 

(Counts/sample  time) 


R 

(Signal/Noise) 


.  100  1 

2.5  5 

10.0  10 

40.0  20 


The  above  tables  show  that  daylight  soundings  of  HF  are  just  feasible 
given  the  scenario  used  in  this  thesis  since  a  signal  count  of  at  least 
36.2  is  needed  in  ten  sample  times.  The  best  that  can  be  done  from 
Table  XI  is  about  36.9  counts.  Nighttime  soundings  look  much  more  pro¬ 
mising.  A  signal  count  of  36  will  yield  a  signal  to  noise  of  about 
six. 


Beam  Divergence 

A  beam  divergence  of  3.33  X  10  ^  radians  was  used  in  this  analysis 
so  that  the  laser  spot  size  would  be  on  the  order  of  the  plume  size 
(about  10  meters).  To  determine  if  the  laser  needs  expanding  optics 
or  not,  it  is  necessary  to  make  an  estimate  of  the  exit  aperature  re¬ 
quired  at  the  laser  to  obtain  the  above  0Q.  Using  equations  (43)  and 
(44),  an  exit  aperature  of  4.32  cm  is  required.  With  this  size  beam, 
expanding  optics  are  probably  necessary. 

Conclusions 

Detection  of  HF  through  LIDAR  techniques  does  appear  to  be  feasible 
if  the  soundings  are  done  at  night.  As  in  the  case  of  HCl,  if  tunable 
lasers  are  developed  that  can  operate  in  the  UV  or  IR  beyond  3  ym,  then 

63 


- nr:*.  _ 


HF  detection  should  be  much  easier.  Sounding  periods  longer  than  1 
second  should  also  contribute  more  signal  counts  to  t lie  LIDAR  receiver 
and  it  is  quite  likely  that  extended  soundings  could  be  carried  out  in 
daylight  with  acceptable  signal  to  noise.  Extended  soundings  were  not 
considered,  as  mentioned  previously,  to  avoid  a  complicated  integration 
over  aspect  angle  (0)  with  time. 
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D.  Diatomic  Iodine  (1^) 

Diatomic  iodine  is  a  molecule  of  interest  because  various  iostopes 
of  molecular  iodine  are  by-products  of  nuclear  explosions.  Thus,  detec¬ 
tion  of  I,  in  a  localized  area  of  the  atmosphere  through  LIDAK  techniques 

could  be  a  promising  way  of  remotely  monitoring  nuclear  testing  activity 

131 

around  the  world.  Ideally,  iodine  molecules  involving  isotopes  I  , 

133  135 

I  ,  or  1  should  be  considered  in  the  LIDAR  analysis,  however,  only 

127 

sufficient  data  on  was  available,  so  this  molecule  was  analysed. 

It  is  believed  that  the  diatomic  iodine  with  various  isotopes  should 
have  comparable  absorption  cross-sections,  but  the  wavelengths  of  the 
various  transitions  will  vary  from  isotope  to  isotope.  The  basic 
scenario  considered  is  the  space  shuttle  LIDAR  platform  in  a  250K11 
circular  orbit  and  the  cloud  lying  in  the  range  of  30,000  to 
50,000  ft.  (9-15KM).  Look  angles  of  0°  (nadir),  30°  and  45°  are  con¬ 
sidered.  A  1  joule  per  pulse,  10  PPS  dye  laser  is  again  considered  and 
a  total  sounding  time  of  one  second  is  assumed  to  allow 0  to  remain 
constant.  A  sample  time  (^t)  of  6.6  psec  is  used.  The  cloud  is 
considered  relatively  hot  at  500°K,  consistant  with  the  expected  con¬ 
ditions  shortly  after  a  nuclear  blast.  The  receiver  aparature  is 
.785  M  (1  meter  in  diameter). 

Transition 

Fortunately,  the  first  few  electronic  transitions  of  Ig  are  in  the 

visible  wavelengths  and  are  thus  accessible  to  the  dye  laser.  The 

3 

electronic  transition  n  is  chosen.  The  Franck-Condon  principal 

favors  the  vibrational  transition  v"=G v'=25  (ref  17).  The  rotational 
fine  structure  transition  of  interest  will  be  determined  by  what  J" 
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level  is  the  most  populated  at  500°K.  Using  equation  (53)  with 

B  =  .03729  cm  '  (rei  13),  we  find  that  J"  =  68.  Therefore,  we  may 

max 

excite  either  (X'X,v"=0,  J"=68 -*■  E3it,v'=25,J  ’=69  or  67)  lines  which  are 
respectively,  the  R(68)  or  P(68)  lines.  The  fluorescence  emission 
lines  of  interest  are: 

B3Tt,v'=25a-X'Z,v"=0 

and: 


1) 

J  '=69->J"=68 

R(68) 

2) 

J'=69>J"=70 

P(70) 

3) 

J'=67+J"=66 

R(66) 

4) 

J'=67-*-J"=68 

P(68) 

The  R(68)  line  is  chosen  arbitrarily  for  laser  excitation.  This  leaves 
a  choice  of  (1)  or  (2)  above  for  the  fluorescence  wavelength.  The  P(70) 
transition  is  chosen  for  fluorescence  since  a  wavelength  other  than 
the  laser  wavelength  is  desired  for  detection  to  eliminate  background 
from  Mie  scattering.  The  laser  and  fluorescence  wavelengths  are  given 
respectively  as: 

Xlaser(cm)  =  (Aex’2->B3)T+AV '=<W-25  +AeJ"6B>J  ,=69) 

and: 

^p(cm)  =  (AeB3u->xi2+Ae:vt_25->v"=0  +  ^ej'=6?+J"=70^ 
where : 

4Ex'^»3,  -4£bV*'I  -  l5769-01  <re£  16) 

and: 

i 

Aev"=0-»u'=25  =Aev’=2W=0  =lev'=25  “ev"-(J 
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where : 
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where  B'  and  B"  represent  respectively  the  rotational  constant  values 
for  the  upper  and  lower  electronic  states.  The.  B  coefficients  were 
determined  from  the  formula: 


B  =  Bc-ae  (v  +  '|-)+Ye^V  +  y)^  +  6e(v  +  '|')^  (68) 

where  the  coefficients  are  correction  factors  to  the  rotational  constants 
due  to  changes  in  the  equilibrium  distance  between  the  two  iodine  atoms 
due  to  variations  in  vibrational  energy.  These  correction  factors, 
like  the  vibration  case,  can  change  depending  on  what  electronic  state 
the  molecule  is  in.  For  the  ground  electronic  state  we  have  (ref  18:40): 

B  "  =  .037389  cm"1 
e 

a  "  =  .000121  cm"1 
e 

The  other  terms  are  neglegible  for  v"=0. 

The  rotational  coefficients  for  the  upper  elctronic  state  are  (ref  17:25): 

B  ’  =  .028969  cm"1 
e. 

a  1  =  .0001562  cm  1 
e. 


v  1  =  4X10~7  cm' 
•e. 


6  '  =  3.5X1 O-8  cm"1 
e 

By  substituting  these  values  into  equations  (67)  and  (68),  we  obtain: 

Agt,i  . 0  T i  ,n  *»  20.046053  cm  1  for  the  laser 

and 


"=59>J’=70 


=  30.423376  cm 


for  the  fluorescence. 


Finally: 


Xlaser  “  05769.01  +  2550.6619  -  20.046053)"1  cm 
Afluor  "  (15769*01  +  2550.6619  -  30.423376)"1  cm 
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o  o 

which  gives  5464. IA  and  5467. 2A  for  the  laser  and  fluorescence  wave¬ 
lengths  respectively.  Since  the  laser  and  fluorescence  wavelengths  are 
very  close  to  each  other,  a  narrower  laser  linewidth  than  that  used 
previously  is  chosen  so  that  there  is  no  chance  of  the  laser  beam 
exciting  other  hyperfine  rotational  lines. 


Fluorescence  Efficiency 


The  fluorescence  efficiency  (S^)  is  the  product  of  the  Franck- 
Condon  factor  and  the  rotational  factor.  From  (ref  19:329)  the  Franck- 


Condon  factor  for  the  V  ,=25'*'\,,a0  transition  is  given  as  (.027).  The 
J  =69  rotational  state  can  decay  into  two  different  lower  state  rota¬ 
tional  levels  J"  =t70  and  J"  =68  with  equal  probability.  Therefore  the 
fraction  of  the  total  number  of  molecules  in  the  v'=25,  J’=69  state  that 
will  decay  to  the  v"=0,  J" =70  level  is: 

0F  =  (.5)  (.027)  =  .0135 


Absorption  Cross-Section 


The  linewidth  of  the  Doppler-broadened  profile  is  given  by 


equation  (24)  with: 

VQ  =  5.4867X1014  HZ 
T  =  500°K 

M  =  2(127  AMU)  ( 1.66X10_24g/AMU)  =  4.2164X10-22  g 


Using  equation  (27)  we  have: 


AvI2  =  5.4965X108HZ 
and 


r  J?L_  °  ft 

=  2.998X10  sec( .  IA)  =  110.46X10  HZ 

(5464.1  A)2 
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Again  we  have  Av^  >:AVj  ,  and  we  may  use  the  approximation  for  the 

effective  cross-section,  equation  (30).  From  (ref  20:1381)  the  life- 

3  1  -7 

time  (x )  of  the  (13  ti-+  X  E  3V  '=25-vv"=0)  transion  is  given  as  7.2X10 

sec.  The  rotational  fine  structure  within  the  vibronic  transitions 
do  not  significantly  perturb  the  A'  value,  and  so  these  small  correc¬ 
tions  are  ignored.  Using  equation  (26)  with: 

8b  2J'+i  2(69)+!  _ 

g  =  2J"+ 1  =  2(70)+ 1  ~ 
a 

and 

o 

X  =  5467.2  A 

we  have: 

°eff  =  1.5252X10-14  cm2  =  1 .5252X1 0_18m2 


LIDAR  Receiver 

The  7265  RCA  photomultiplier  tube  seems  to  be  a  good  choice  for 

detection  of  the  fluorescence  (ref  5:211).  The  photocathode 

(Na  Sb  K  Cs)  has  a  relatively  high  responsivity  and  a  quantum  efficiency 
o 

at  5461  A  of  about  (.1).  The  throughput  factor  (e)  is  the  same  as 

that  used  for  sodium.  The  reflectivities  of  the  mirrors  within  the 

o 

detector  telescope  have  basically  the  same  values  at  5461  A  as  they  had 

o 

for  the  sodium  wavelength  at  5896  A. 


Transmission 

A  table  of  transmissivity  vs.  penetration  depth  (\)  into  the 

atmosphere  was  generated  by  interpolation  between  transmission  values 
o  o 

for  the  5461  A  and  6328  A  laser  lines  in  the  same  manner  as  explained 
in  the  sodium  section.  Finally,  a  table  of  transmission  vs.  penetration 
depth  into  the  atmosphere  was  generated  (Table  XIV).  As  in  the  case 
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of  sodium,  penetration  depth  (x)  is  defined  as  zero  kilometers  at  100  KM 


above  the  earth  and  100  KM  at  ground  level. 

Table  XIV 

Transmission  vs.  Penetration 


o 


Depth  at 

5462  A 

X 

T, 

X 

T_ 

X 

T_ 

30 

.9992 

80 

.9896 

91 

.9352 

50 

.9990 

81  • 

.9868 

92 

.9275 

55 

.9939 

82 

.9834 

93 

.9193 

60 

.9987 

83 

.9796 

94 

.9106 

65 

.9983 

84 

.9755 

95 

.9007 

70 

.9973 

85 

.9711 

96 

.8864 

75 

.9967 

86 

.9662 

97 

.8590 

76 

.9958 

87 

.9609 

98 

.7991 

77 

.9947 

88 

.9551 

99 

.6654 

78 

.9934 

89 

.9489 

100 

.4096 

79 

.9917 

90 

.9423 

It  is  assumed  that  the  natural  concentration  of  l„  in  the  atmosphere  is 
small  enough  so  that  losses  due  to  resonance  absorption  (K^)  is  nearly 
zero. 


Noise  Analysis 

The  noise  analysis  follows  the  same  general  procedure  as  outlined 

in  the  sodium  case.  Once  again  the  curves  in  ref  4:42  were  used  to 

determine  radiances  of  reflected  sunlight  and  moonlight  off  of  the  earth 
o 

at  5461  A.  This  fluorescence  wavelength  does  not  lie  on  a  Fraunhofer 
absorption  line;  thus  noise  values  for  I ^  will  be  greater  than  those 
for  sodium.  The  radiances  of  reflected  sunlight  and  moonlight  at  5461 
are  respectively: 

17  0  2  ‘ 

L£  =  1.7424X10  photons/scc  A  str  m 

and 

12  0  2 
“  2.1818  X10  photons/sec  A  str  m 
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Returning  Co  equation  (40)  with: 

A  =  .785  m2 

0^  =  10  ^  rad 
u  c 

AA  ~  -5  A 
n  =  .i 
e  =  .765 
At=  6.6  psec 

The  number  of  noise  photons  detected  per  sample  time  is: 

Ng  =  278.29  photons/sample  time 
-3 

Nm  =  3.48X10  photons/sample  time 


The  noise  due  to  dark  current  can  again  be  neglected  if  the  photocathode 
is  cooled  to  220°K.  Using  equation  (41),  one  can  determine  the  minimum 
number  of  signal  photons  (Nmj,n)  required  to  give  a  particular  signal  to 
noise  ratio  (R)  for  daylight  and  night  observations.  Using  the  standard 
LIDAR  values  of  laser  pulse  frequency  (fp)  equal  to  10  Hertz,  LIDAR 
sounding  time  (x^)  of  1  second  and  dynode  gain  (6)  of  3.38,  Tables  XV 
and  XVI  present  Nm^n  vs.  R  first  for  daylight  observations  and  then 
nighttime . 

Table  XV 

Iodine  N  .  vs.  R 
mm 


(Daylight) 


N  . 
mm 


6.36 

12.86 

26.31 

70.37 

157.31 

207.70 

Note  that  N  .  (I,)>N  ■  (Na), 

mm  2  mm  ’ 


R 

1 

2 

4 

10 

20 

25 

this  is  mainly  due  to  the  Fraunhofer 
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attenuation  effect  which  is  present  for  sodium  but  not  for. iodine. 


Table  XVI 

Iodine  N  .  vs.  R. 
rmn 

(Nighttime) 


» 

min 

R 

.  145 

1 

.572 

2 

2.28 

4 

14.2 

10 

56.83 

20 

88.79 

25 

Note  that  the  nighttime  values  are  comparible  to  those  of  sodium.  This 
is  because  the  noise  terms  at  night  are  so  small  that  they  can  be 
neglected.  The  only  appreciable  noise  is  the  signal  noise  and  this  is 
the  same  for  both  and  sodium.  By  choosing  a  desired  signal-to-noise 
value  (R),  one  can  express  the  minimum  column  density  required  as  a 
function  of  range,  look  angle  and  laser  energy  per  pulse  by  using 
equation  (1).  For  R=10,  equation  (1)  yields  “  70.37.  Tables 

of  minimum  column  density  (Km^n)  vs.  laser  energy  per  pulse  can  be 
constructed  for  daylight  and  night  soundings  (Tables  XVII  and  XVIII). 

Table  XVII 

Daylight  Minimum  Detectable  Column  Density 
vs.  Laser  Energy  -  Iodine 


(Joules/pulse) 


(10 


mm 

-16  -2, 
m  , 


1.0 

1.52 

.8 

1.90 

.6 

2.54 

.4 

3.81 

.2 

7.63 

.  1 

15.2 

.01 

152.0 

.002 

762.0 
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Table  XVIII 


Night  Minimum  Detectable  Column 
Density  vs.  Laser  Energy  -  Iodine 


I 

Joules/pulse 


K  . 
mm 


) 


1.0 
.8 
.6 
.4 
.2 
.  1 
.01 
.002 


.307 

.385 

.511 

.770 

1.54 

3.07 

30.7 

154.0 


Tables  XVII  and  XVIII  utilized  the  following  values: 


At  =  6,6  usee  0  =  .0135 

r 

AZ  =  1KM  R  =  10 


r)£  =  .0765  Shuttle  altitude  =  250  KM 

Z  -  235  KM  a  -  1.52X1018  m2 

0  =  0° 

A  table  of  minimum  column  density  vs.  range  and  look  angle  for  daylight 
and  night  observations  is  presented  next  in  Table  XIX.  The  parameters 


not  varied  are  the  same  as  the  ones  used  in  Tables  XVII  and  XVIII.  The 


laser  energy  per  pulse  was  fixed  at  one  joule. 


Table  XIX 


Minimum  Detectable  Column  Density  vs.  Range 
and  Look  Angle  -  Iodine 


z , 

Z 

e 

(T  T  feC  ° 

Say 

Sight 

V 

(KM) 

(KM) 

(Deg) 

r  t 

( 10 16m"2) 

( )015>f  2) 

150 

100 

0 

1.0 

.26 

.53 

175 

175 

0 

.998 

.80 

1.61 

200 

200 

0 

.998 

1.04 

2. 10 

215 

215 

0 

.996 

1.21 

2.44 

235 

235 

0 

.943 

1.52 

3.07 

237 

237 

0 

.924 

1.58 

3.  19 

240 

240 

0 

.887 

1.69 

3.41 

243 

243 

0 

.846 

1.81 

3.66 

245 

245 

0 

.812 

1.92 

3.89 

248 

248 

0 

.640 

2.50 

5.04 

150 

173.21 

30 

1.0 

.78 

1.58 

175 

202.07 

30 

.998 

1.06 

2. 14 

200 

230.94 

30 

.998 

1.39 

2.80 

215 

248.26 

30 

.995 

1.61 

3.25 

235 

271.35 

30 

.934 

2.05 

4.15 

237 

273.66 

30 

.913 

2.13 

4.30 

240 

277.13 

30 

.871 

2.29 

4.63 

243 

280.59 

30 

.824 

2.49 

5.00 

245  . 

282.90 

30 

.786 

2.65 

5.33 

248 

286.37 

30 

.597 

3.57 

7.22 

150 

212,13 

45 

1.0 

1.17 

2.36 

175 

247.49 

45 

.997 

1.60 

3-22 

200 

282.84 

45 

.997 

2.09 

4.22 

215 

304.06 

45 

.994 

2.42 

4.89 

235 

332.34 

45 

.920 

3.12 

6.30 

237 

335. 17 

45 

.894 

3.27 

6.59 

240 

339.41 

45 

.844 

3.55 

7.  15 

243 

343.65 

45 

.789 

3.89 

7.85 

245 

346.48 

45 

.745 

4.19 

8.44 

248 

350.72 

45 

.532 

6.00 

12.1 
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Note  that  it  is  assumed  that  T  ='J’  which  is  a  good  approximation  since 
the  wavelengths  differ  by  only  3  Angstroms.  Finer  divisions  of  Z  between 
235-243  KK  were  made  since  the  cloud  is  considered  to  lie  at  an  altitude 
of  9-15  Km  above  the  earth. 

Beam  Divergence 

-4 

A  laser  beam  divergence  of  10  radians  was  considered  adequate  for 

L1DAR  sounding  of  (ref  4:31).  To  confirm  this,  one  can  solve  equation 

(35)  with  Z  equal  to  350  KM  which  is  the  worse  case  considered  here. 

This  yields  a  spot  size  diameter  of  35  meters  which  is  probably  smaller 

-4 

than  the  extent  of  the  cloud.  With  Q^equal  to  10  radians,  one  can 
calculate  the  required  exit  aperature  of  the  LIDAR  transmitter  using 
equation  (36).  The  aperture  diameter  turns  out  to  be  6.7  millimeters. 
Therfore,  no  expanding  optics  are  necessary  because  the  beam  exiting  the 
laser  will  probably  be  just  a  few  millimeters  v?ide. 

Conclusions 

In  general,  the  required  column  densities  are  a  thousand  times  great¬ 
er  for  diamotic  iodine  than  they  are  for  sodium.  However,  they  are  still 
very  small  densities  and  detection  of  1^  after  a  nuclear  explosion  seems 
to  be  promising. 


E.  Diatomic  Hydrogen  (l^) 

Diatomic  hydrogen  is  of  interest  because  its  detection  around  an 
industrial  complex  may  indicate  possible  activities  involving  nuclear 
processes.  Unfortunately  again,  electronic  transitions  of  are  too 
far  out  in  the  ultraviolet  to  be  accessible  to  the  Dye  laser.  Vibration- 
rotation  levels  are  the  only  possibilities  for  dye  laser  excitation. 
Diatomic  hydrogen  has  no  permanent  dipole  moment,  and  so  only  quad- 
rapole  transitions  are  allowed.  The  lifetimes  of  these  transitions  are 
on  the  order  of  10^  seconds  (ref  22:133)  which  means  that  the  effective 
absorption  cross-section  is  roughly  six  orders  of  magnitude  smaller 
than  that  of  HCl.  Since  LIDAR  detection  of  HC1  is  marginal  at  best 
because  of  it's  absorption  cross-section,  it  is  clear  that  LIDAR  detection 
of  1^  through  vibrational-rotational  excitation  is  impossible  unless  a 
tunable  ultraviolet  laser  is  developed  to  access  the  larger  absorption 
cross-sections  of  electronic  transitions. 

F.  Hydrogen-Tritium  (HT) 

HT,  like  H^,  is  of  interest  because  it  may  also  indicate  nuclear 

activity  if  detected  in  an  industrial  area.  Unfortunately,  LIDAR 

detection  of  this  molecule  seems  unlikely  because  it  probably  has  a 

very  small  permanent  dipole  moment.  This  conclusion  was  reached  after 

discovering  that  HD  (hydrogen-deuterium)  has  a  permanent  dipole  moment 

four  orders  of  magnitude  smaller  than  HCl  (ref  16:259).  Since  HT  and 

HD  differ  by  only  one  neutron,  it  is  suspected  that  they  have  similar 

dipole  moments.  If  this  is  true,  then  in  all  probability  the  dipole 

moment  transition  matrix  elements  (|RV  j)  are  very  small  since  these 

V",J" 

matrix  elements  depend  on  M(r)  (see  equation  57).  Since  the  absorption 
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cross-sect  ion  is  proport ional  to  the  square  of  the  transition  matrix 
element,  it  too  will  be  small.  Therefore,  it  appears  that  detection 
of  HT  is  impossible  with  present  day  tcchnoJogy.  Electronic  transitions 
again  are  too  far  out  in  the  ultraviolet  for  today's  tunable  lasers. 

G.  Xenon  (Xe)  and  Krypton  (Kr) 

Xenon  and  krypton  are  possible  by-products  of  a  nuclear  blast. 
Therefore,  they  are  of  interest  for  the  same  reason  as  diatomic  iodine. 
Unfortunately,  these  elements  cannot  be  detected  with  a  tunable  dye 
laser  because  the  transition  from  the  ground  state  to  the  first  excited 
state  for  both  atoms  is  in  the  ultraviolet.  Of  course,  since  krypton 
and  xenon  are  atoms,  there  is  no  vibrational -rotational  spectra  to 
excite.  As  mentioned  before,  a  tunable  dye  laser  that  can  operate  in 
the  ultraviolet  is  the  only  means  of  conducting  LIDAR  resonance  scatter¬ 
ing  techniques  on  Xenon  and  krypton. 

H.  Boron  Oxide  (B2®3^’  Bor^c  Acid  (HBC^)  and  Aluminum  Oxide  (A^O^) 

These  molecules  are  interesting  because  their  detection  in  a  rocket 
plume  may  allow  identification  of  the  rocket  fuel  being  used  and  per¬ 
haps  allow  determination  of  the  vehicle  performance.  Very  little  data 
was  found  on  these  molecules  especially  in  the  absorption  cross-section 
area,  which,  from  the  work  of  this  thesis,  appears  to  be  the  most  import¬ 
ant  parameter  to  determine.  Some  data  was  gathered  on  the  molecular 
spectra  of  gaseous  HRC^  and  B2®3  an<^  data  is  presented  in  Tables 

XX  and  XXI  for  possible  use  for  anyone  following  up  on  this  study 
(ref  23:269)  and  (ref  24:298).  No  significant  data  on  aluminum  oxide 
was  found. 
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Table  XX 


\ 

0 


Band 

HBO2  Spectral  Data 

Band  Center 

(cm  ') 

B=0  Stretch 

4.97 

pm 

0-H  Stretch 

2.27 

pm 

B 

Table  XXI 

2^3  Spectral  Data 

Spectral  Location  (cm 

Band 

B2“03 

Vj(B=0  Symmetric  Stretch) 

2087 

2021 

V„(B=0  Symmetric  Stretch) 

760 

757 

V^(0-B-0  Symmetric  Bend) 

498 

478 

v^(B-0-B  Bend) 

236 

235 

V^(B=0  Asymmetric  Stretch) 

2095 

2029 

V^(B-0  Asymmetric  Stretch) 

1270 

1268 

v  (O-B-0  Asymmetric  Bend) 

8 

568 

519 

Vg(Out  of  Plane  Bend) 

493 

478 

7  9 


I.  Tritiated  Water  (HTO) 

Like  and  HT,  detection  of  tritiated  water  around  an  industrial 

complex  may  indicate  nuclear  activity  in  the  area.  Very  little  data 

was  obtained  on  HTO,  and  so  an  extended  analysis  of  this  molecule  was 

not  possible.  However,  it  may  be  possible  to  use  absorption  cross- 

sectional  data  of  ordinary  H^O  and  utilize  this  for  HTO.  From  (ref  6:106) 

o 

absorption  cross-section  values  in  the  visible  (6944  A)  range  of  the 

.  -23  -2 

spectrum  were  found  to  lie  m  the  range  of  10  cm  .  This  is  not  a 
particularly  impressive  cross-section  and  it  appears  that  HTO  sounding 
in  the  visible  spectral  regime  will  be  very  difficult.  Spectral  data 
was  obtained  (ref  25:454)  and  the  findings  are  tabulated  in  Table  XXII. 

Table  XXII 


HTO  Spectral  Data 


Band 

Bending 

Symmetric  Stretch 
Asymmetric  Stretch 
Asymmetric  Stretchf 
Symmetric  Stretchf 
Bendingt 
Bendingtt 
Sym-Asym  Stretch 
Bendt-Asym  Stretch 
Sym  Stretch-Bendf 
Bend-Asym  Stretch 

t  -  1st  overtone 


Spectral  Location  fcm  *) 

1324 

2300 

3711 

7221 

4565 

2643 

3933 

6022 

6322 

4929 

5032 


ft  -  2nd  overtone 


\  V.  Cone lus i ons  and  Recommendations 

V  . . -  ■  ■ - ■- — - -  ■  ■-  — 

It  appears  that  L1DAR  detection  o£  sodium, hydrogen  fluoride  and 
diatomic  iodine,  is  very  promising.  This  is  due  primarily  to  the  large 
absorption  cross-sections  that  the  three  species  possess.  The  absorp¬ 
tion  cross-section  was  found  to  be  the  main  parameter  determining  if  a 
species  could  or  could  not  be  detected  through  LIDAR  techniques. 

Diatomic  hydrogen  and  hydrogen  tritium  had  very  small  vibrational- 
rotational  absorption  cross-sections  and  thus  were  unlikely  candidates 
for  LIDAR  sounding.  Hydrogen  chloride  had  a  better  cross-section  but 
was  still  small  enough  to  say  that  its  detection  through  resonance 
fluorescence  techniques  would  be  marginal  at  best.  Xenon  and  krypton 
are,  of  course,  not  molecules;  thus  they  had  no  vibration-rocation 
spectra.  They  must  be  excited  electronically.  The  transition  from 
ground  state  to  the  first  excited  electronic  state  unfortunately  lies 
too  far  out  in  the  ultraviolet  and  beyond  the  capabilities  of  a  dye 
laser.  They  therefore  cannot  presently  be  detected  through  the 
resonance  technique. 

Cross-sectional  data  on  boric  acid,  boron  oxide  and  aluminum  oxide 
could  not  be  found,  and  so  these  molecules  were  not  analysed  in  detail. 
Finding  these  cross-sections  are  essential  if  future  work  is  to  be  done 
in  this  area.  The  development  of  a  tunable  ultraviolet  laser  would 
make  detection  of  most  of  these  atoms  and  molecules  feasible  because 
one  could  then  access  the  relatively  hugh  electronic  absorption  cross- 
sections  of  the  particles.  Einstein  (A1)  coefficients  for  hot  band 
transitions  (Av=l,v"=0)  may  be  useful  information  to  obtain  because 
some  of  these  transitions  have  higher  absorption  cross-sections. 
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Refinements  of  work  presented  here  would  involve  using  computer 
codes  such  as  HITRAN  developed  by  AFGL  to  obtain  much  more  precise 
transmission  coefficients  for  the  laser  and  fluorescence  lines.  Longer 
sounding  scenarios  should  be  examined  to  see  if  signal  to  noise  can  be 
improved  significantly.  This  will  require  a  small  computer  code  to 
account  for  the  changes  in  look  angle  and  range . 
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Hemote  sensing  from  the  space  shuttle  utilizing  laser-induced  fluorescence  tech¬ 
niques  was  investigated.  A  feasibility  study  of  LIDAR  (Light  Detection  and 
Ranging)  techniques  to  detect  twelve  gases  of  scientific  and  military  interest 
was  carried  out.  Evaluation  of  a  European  Space  Agency  LDIAR  hardware  design 
was  also  a  task.  The  main  emphasis  was  on  determination'  of  background  noise, 
spectral  transitions,  atmospheric  transmission  and  molecular  abqsi’ption  cross- 
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required  to  be  detectable.  Atomic  sodium,  molecular  iodine  and  hydrogen 
flouride  were  determined  to  be  promising  candidates  for  LIDAR  resonance 
detection.  Other  molecules  investigated,  such  as  hydrogen  chloride,  xenon 
and  krypton  were  considered  marginally  detectable  at  best  because  of  their 
small  absorption  cross-sections.  Lack  of  sufficient  data  prevented  extended 
analysis  of  the  remaining  molecules. 
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